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ABSTRACT 


Calculations  are  described  which  yield 
absorption  line  parameters  for  the  three  iso- 
topes of  molecular  oxygen:  leO?1  lc0iR0,  and 
1r02,  in  the  format  of  the  AFCRL  Atmospheric 
Absorption  Line  Parameters  Compilation.  The 
line  parameters  are:  transition  frequency, 
integrated  line  strength  at  296K,  line  width, 
lower-state  energy,  and  identifying  quantum 
numbers.  These  parameters  are  required  as 
input  to  the  SLAM  program  described  elsewhere 
which  calculates  the  attenuation  vs.  altitude 
at  any  fixed  frequency  in  the  mi llimeter-to- 
aubnullimeter  wave  region. 

The  calculations  were  performed  on  the 
XDS  Sigma  9 computer  at  RRI  by  means  of  APL 
programs  which  ate  listed  and  explained  in  the 
report.  Several  sets  of  published  molecular 
constants  were  used  to  obtain  an  estimate  of 
the  degree  of  reliability  of  the  resulting  line 
parameters.  Appendices  include  the  predicted 
transition  frequencies,  transition  moments,  and 
lower  state  energies  for  each  set  of  molecular 
constants  employed;  tabulations  of  rotational 
state  sums  vs.  temperature;  line  parameters  for 
each  of  the  three  isotopes  studied,  separately; 
and  a list  of  the  318  absorption  lines  below 
300  cm  1 whose  strengths  exceed  3.7  E-30  cm-1 
per  molecule  cm"2  at  296K. 
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THE  CALCULATION  OF  MILLIMETER  AND  SUBMILLIMETER 
WAVE  ABSORPTION  LINE  PARAMETERS  FOR  THE  MOLECULAR 
OXYGEN  ISOTOPES:  l60,  . 160la0.  AND  leO. 


Introduction 

This  Technical  Report  describes  a portion  of  the  calculations 
performed  at  RRI  on  the  atmospheric  propagation  characteristics 
of  electromagnetic  radiation  of  millimeter  and  submillimeter 
wavelengths.  The  calculations  described  here  are  those  of  a 
set  of  absorption  line  parameters  for  three  isotopes  of  molecu- 
lar oxygen:  160^,  leoi80,  and  le02.  These  line  parameters  are 
required  by  the  SLAM  program  (described  in  a separate  Tech.  Re- 
port) to  calculate  the  attenuation  vs.  altitude  at  any  fixed 
frequency  in  the  covered  wavelength  region.  The  line  parame- 
ters are:  frequency  (cm  i),  strength  at  296k  (cm  1 per  molecule 
per  cm'),  linewidth  (cm  1 per  atm),  and  lower  state  energy 
(cm-1)  for  each  transition,  labelled  by  the  appropriate  quantum 
numbers,  date  of  calculation  (month,  year),  isotope,  and  species 
(7  for  oxygen).  The  calculations  were  done  on  the  basis  of 
several  sets  of  published  molecular  constants1  7 by  means  of 
AFL  programs  written  for  the  XDS  Sigma  9 computer  at  RRI.  The 
programs  are  listed  and  explained  in  this  report. 

The  succeeding  sections  of  this  report  discuss  the  molecular 
theory  of  the  oxygen  molecule,  the  calculation  of  line  positions, 
lower  state  energies,  and  integrated  line  strengths,  evaluation 
of  line  widths,  and  the  format  of  the  line  parameters  listed  in 
the  Appendices.  In  addition  to  the  line  parameter  listings  (for 
ail  isotopes  except  i80i70),  the  Appendices  contain  the  following: 
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a brief  summary  of  APL  notational  conventions,  predictions 
of  transition  frequencies  and  lower  state  energies  to 
accuracy  (both  in  GHz  and  in  cm”1)  for  each  set  of  molecular 
parameters  employed,  a tabulation  of  rotational  state  sums  vs. 
temperature,  and  several  other  listings  of  general  interest, 
but  not  required  in  the  body  of  the  report. 

Molecular  Theory  of  the  Oxygen  Isotopes 

The  first  essentially  correct  theoretical  treatment  of  the 
fine  structure  of  the  molecular  oxygen  ground  state  was  that  of 
Tinkham  and  Strandberg.1  This  theory  (applicable  to  1602  and 
1802)  predicted  that,  in  addition  to  the  well-known  microwave 
spectrum,  oxygen  also  possessed  a "forbidden"  submillimeter  spec- 
trum, and  made  predictions  of  the  line  strengths.  The  essential 
correctness  of  the  theory  (with  the  exception  of  the  strength 
prediction  for  one  type  of  "forbidden"  transition)  was  shown 
by  Gebbie,  et  al.,8  on  the  basis  of  laboratory  measurements  at 
low  resolution,  as  well  as  by  the  observation  of  atmospheric 
absorption  at  high  altitudes  using  the  Sun  as  a source.91 10 
In  his  recent  thesis,  Steinbach  has  succeeded  in  clarifying  the 
theory  with  the  aid  of  the  transformation  theory  for  spherical 
tensor  operators, 11  and  it  is  based  on  Steinbach ‘s  formulation 
that  the  RRI  calculations  of  line  strengths  have  been  performed, 
especially  since  it  places  le0lB0  on  an  equal  footing  with  the 
homonuolear  isotopes. 7l  11 

The  basic  physical  features  are  as  follows:  None  of  the 
oxygen  molecular  isotopes  possess  a permanent  electric  dipole 
moment.  For  the  le02  and  l802  isotopes,  this  follows  from  the 
homonuclear  symmetry,  whereas  for  l60l80  (as  well  as  le0l70), 
this  is  an  experimental  fact.11  The  electronic  ground  state 
being  3Zg  , the  average  (electronic)  orb j cal  angular  momentum 
vanishes  (A  — 0)-— although  there  is  an  instantaneous  non— zero 
value  of  the  orbital  angular  momentum  which  produces  a precessing 
magnetic  dipole  moment  of  orbital  origin.  Since  the  electronic 


-2- 


T -1/306-3-14 


RIVERSIDE  RESEARCH  INSTITUTE 


ground  state  is  a triplet  state,  the  total  electronic  spin  is 
S = 1,  so  that  molecular  oxygen  has  a permanent  magnetic  dipole 
moment  of  approximately  two  Bohr  magnetons.1  The  oxygen  mole- 
cule interacts  via  this  dipole  moment  with  the  magnetic  field 
associated  with  an  electromagnetic  wave.  It  has  been  noted  by 
several  authors1* 1 1 * 1 that  oxygen  is  capable  of  very  strong 
atmospheric  absorption  of  electromagnetic  radiation  despite  the 
weakness  of  magnetic  dipole  transitions  relative  to  electric 
dipole  transitions  (,by  two  or  more  orders  of  magnitude,  in  gene- 
ral) by  virtue  of  its  abundance  (21^  by  volume  of  air  at  all 
altitudes).  (By  way  of  contrast,  water,  with  the  large  elec- 
tric dipole  moment  of  1.9  Debye  units,  decreases  rapidly  in 
concentration  with  altitude  above  3 ) 

As  is  typical  for  microwave  and  submillimeter  wave  absorp- 
tion, the  molecular  oxygen  absorption  frequencies  in  this  region 
are  associated  with  changes  in  the  rotational  energy  of  the  mole- 
cule, with  no  change  taking  place  either  of  the  vibrational  or  of 
the  electronic  state.  (At  room  temperature,  both  the  ground 
vibrational  state  and  the  first  excited  vibrational  state  are 
significantly  populated  and  are  taken  into  account  in  these 
calculations.)  The  energy  in  any  given  rotational  state  is  the 
sum  of  three  terms:  the  pure-rotation  term,  HrQt  = BN2  (where 
N is  the  rotational  angular  momentum  operator  describing  end- 
over— end  rotations  of  the  molecule,  and  B is  the  rotational 
parameter — itself  a function  of  N when  stretching  , i.  e., 
rotation-vibration  interaction,  is  taken  into  account);  a spin- 
rotation  (or  spin-orbit1)  term,  Hg_r  = u.N*S  (where  S is  the 
total  electron  spin  operator,  and  n.  is  the  spin-rotation  coupl- 
ing parameter — likewise  N-dependent  under  centrifugal  stretch- 
ing) ; and  a spin-spin  term,  Hs_s  = ( 2/3)  A(  3SZ‘ -S?)  (where  Sz  is 
the  body- fixed  component  of  S along  the  internuclear  axis  of 
the  molecule,  and  X is  the  spin-spin  interaction  parameter,  com- 
parable in  magnitude  to  B) . A spectroscopic  convention  is  fol- 
lowed here:  not  only  frequencies,  but  also  energies  and  rotational. 
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spin-rotation  and  spin-spin  interaction  parameters,  are  expressed 
in  frequency  units:  GHz  or  cr.i-1.  Note  that  the  speed  of  light 
is  29.9792458(1.2)  GHz  per  cm-1,13  so  that  1 cm-1  as  30  GHz. 

Since  the  nuclear  spins  of  l60  and  of  l80  are  zero  (I  = 0), 
nc  hyper  fine  energy  term  appears  in  the  Hamiltonian.  For  “■  O, 
however,  I = 5/2,  enormously  complicating  the  spectrum  of  00 
relative  to  that  of  l60l80,  le02,  or  le02. 14-16  The  hyperfine 
spectrum  of  l60l70  (as  well  as  the  Zeeman  effect  in  the  presence 
of  a d-c  magnetic  field,  for  all  isotopes11)  will  be  treated  in 
another  Technical  Report, 17  based  on  some  results  by  Steinbach. 16 
The  relative  abundances  of  the  various  isotopic  species  of 
molecular  oxygen  are  given  in  Table  I. 

Therefore,  restricting  our  attention  to  oxygen  molecules 
with  1=0,  there  are  only  three  angular  momentum  operators  to 

consider:  N,  the  rotational  angular  momentum,  S,  the  spin  angu- 

^ 

lar  momentum,  and  J = N + S,  the  total  angular  momentum.  J 
enters  because  the  spin-rotation  term  may  be  rewritten  as  a 
function  of  N2,  Is2,  and  J2  alone,  since  N*S  = (l/2)  (N2+S2-J2) . 

As  is  well-known, lf 8* l2,  18  the  ordinary  microwave  absorption 
spectrum  of  le02  can  be  explained  to  a good  first  approximation 
by  assuming  the  spin  tightly  coupled  to  the  rotation  in  a pure 
Hund's  case  (b)  fashion.  However,  such  a coupling  scheme  fails 
to  predict  the  existence  of  the  submillimeter  absorption  lines 
and  also  fails  to  reproduce  all  the  details  of  the  observed 
microwave  spectrum. lf 19  Consequently,  Tinkham  and  Strandberg 
performed  the  energy  calculations  in  a Hund's  case  (a)  basis 
(in  which  J,  N,  S,  and  Mg  are  good  quantum  numbers)  and  trans- 
formed to  a Hund's  case  (b)  basis  (in  which  J,  MJ#  N,  and  S are 
good  quantum  numbers).1  The  calculations  are  fairly  involved, 
and  the  expressions  needed  to  calculate  transition  strengths 
contain  typographical  errors,  but  (with  the  exception  noted 
by  Gebbie  et  al, 8 and  independently  by  Steinbach11)  their  numeri- 
cal values  for  the  transition  matrix  elements  (line  strengths) 
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Table  I:  Relative  Abundances  of  the  Isotopes  of  Op 
Isotopic  Species  Relative  Abundancea 


le0l60 

0.99519 

l60le0 

4.07.10"3 

160170 

7.38-10-4 

18q18q 

4.16-10"6 

I70ie0 

1.51*10"° 

17q170 

1.37-10"7 

Based  upon  the  following  isotopic  abundances  of  atomic  oxygen 

l60:  99-759#,  l70:  0.037#,  l80:  0.204# 

Reference:  Handbook  of  Chemistry  and  Physics,  52nd  Edition, 

1971-1972,  Chemical  Rubber  Co.,  Cleveland,  Ohio. 
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agree  with  those  of  later  workers.11  A different  approach  has 

been  taken  by  Steinbach7* 11  based  on  earlier  (but  incorrect1*') 

calculations  by  Hill  and  Mizushima. 19  In  this  approach,  all 

calculations  are  carried  out  in  a Hund's  case  (b)  basis,  and 

one  finds  the  matrix  elements  of  S in  this  basis.  Steinbach11 

z 

performed  the  calculations  by  using  the  theory  of  irreducible 
tensor  operators  as  summarized  by  Ccok  and  DeLucia13  to  obtain 
explicit  algebraic  expressions  for  the  eigenvectors  describing 
the  stationary  energy  states: 

U'>  = °l  J,  Mj,  N~  J-lj  S>  + 0\  J,Mj,N=J,  S>  +y\  J,Mj,N=J+l,S>  , 

...(1) 

in  terms  of  the  three  parameters,  a,  0,  and  7,  using  the  Wigner 
6-j  symbols.  The  fact  that  (with  one  exception)  either  a or  7 
is  neatly  1 while  the  other  is  almost  (but  not  quite)  zero11 
means  that  oxygen  is  "almost"  Hund's  case  (b),  but  that  N has 
ceased  to  be  a "good  quantum  number"  except  for  "diagonal"  states 
(corresponding  to  aD  = 70  = 0,  0o  = 1 in  Steinbach 's  notation11), 
as  noted  earlier  by  Mizushima  and  colleagues3*4*6*19  and  by 
Tinkham  and  Strandberg.1  In  these  calculations,  the  molecular 
rotational,  spin-spin,  and  spin-rotation  parameters,  B,  A,  and 
M-i  are  assumed  to  depend  as  follows  on  N and  J:3*6*11*19 

B = B0  + BxN(N+1)  + B2N2( N+l)  3 + ••• 

A = Ao  + AiN(N+l)  + ••• 

M.  = M-o  + p.iN(N+l)  + •** 

A = A0  + Ai( J2+J+l)  + •••  for  N f N' , N = J+l 

matrix  elements. 

...(2) 

The  seven  parameters,  B0,  Bi,  B2,  Ao,  Ax,  \i0,  and  Mi  are 
obtained  by  fitting  observed  transition  frequencies  to  the  model 
Hamiltonian,  supplemented  by  calculations  of  the  isotopic 


_ for  N = N' 

matrix  elements. 
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dependence  of  these  parameters,  where  required, 7'  11' 21  The 
separation  between  the  lowest  rotational  levels  in  the  ground 
vibrational  state  (v  = 0)  and  in  the  first  excited  vibrational 
state  (v  = l),  denoted  by  AG1/2  by  Herzberg, l8' 21  is  taken  to  be 
1556.378(7)  cm  1,  based  upon  the  re-evaluation  by  Albritton,  et 
al5  of  electron  optical21  and  microwave3'4  data. 

These  parameters,  as  given  by  Refs.  1 through  7,  have  beer 
inserted  into  APL  "functions"  (programs)  such  as  PARAMSTEIN, 
listed  in  Pig.  1.  A brief  description  of  the  APL  programming 
language  is  contained  in  Appendix  A.  In  Appendix  B,  the  para- 
meters are  explicitly  given,  together  with  the  transition  matrix 
elements,  transition  frequencies,  etc.,  calculated  from  them  as 
discussed  in  the  ensuing  sections.  Listings  of  two  similar  APL 
parameter-setting  programs  used  in  the  RRI  calculations,  PARAM- 
STEINUPPER  and  PARAMETERS,  are  given  in  Appendix  C, 

The  parameters  quoted  for  Tinkham  and  Strandberg1  are  based 
on  the  correspondences:  B(q)t_s  - B0,  4<2B(o)t_s  = Bi,  B2  = 0, 

X(0)T-S  = Xo'  ^*2AlT-S  = Al>  ^0  T-S  = M-o,  M-i  = °*  (See  Ref.  1, 
footnote  4l.)  The  parameters  quoted  for  Albritton,  et  al5  are 
based  on  electronic  optical  data21  for  v = 0,  and  on  ground-state 
microwave  data3' 4 and  upper-state  electronic  optical  data21  for 
v = 1.  (See  Ref.  5,  p.  il6,  Table  IX.) 

Line  Positions  and  Lower  State  Energies 

Except  for  an  arbitrary  additive  constant,  the  eigenvalues 
of  the  Hamiltonian  operator  in  the  stationary  states  | Tp^)>  are 
the  allowed  energy  levels.  (Different  authors1'3'5  choose  dif- 
ferent additive  constants.)  We  follow  the  calculational  scheme 
of  Mizushima,  et  al,3*4'6  and  Steinbach  and  Gordy7* 11  in  which 
the  eigenvalues  are  obtained  by  solving  the  secular  determinan- 
tal  equations: 

I H11  - E I = 0,  H11  = < J=0,Mj,N=l,S  I H | 0,Mj,l,S  >,  ...(3a) 


t-i/306-3-14 
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V PAHAMSTEIN  -.VECBOiVECBl ; VECB  2;  VECLAh  0 ; VECLAM1  ; VECMU  0 ; VECMU1 
C 1 J * 'REFNO ' IS  7,  CORRESPONDING  TO  h.  STEIN BACH  AND  W.  GORDY,  PHYS. 

[2]  ^ REV.  A.  VOL.  11,  NO.  3,  PP . 72  9 TO  731,  MARCH  1 975  . 

[3]  VECB 0*4 3. 10046  40.707409  38.31373 

[4]  VECB1 -~0 . 00014  501  "0.000129  "C. 0001  1 5 

[5]  VECB 2-0  0 0 

[6]  VECLA  .V0<- 59. 501341  59.499097  59.4  96690 

[7]  VECLAM1-  5 . 84  aff”  5 5.312P"5  5.211P"S 

[8]  VECMU  0-~ 0 .2  52  586  "0.2384  08  "0.22443  9 

[9]  VECMUl-~2  .47P"7  "6 . 1 9P"7  "3.51P"? 

[10]  REFNO- 7 

[11]  IS 0-0  0 60  0 0 

[12]  ISOTEXT- 3 7o  ’ ' 

[13  ] isorsmi ; ]*  ’oi  0=010 ' j/sorrm  2;  ]-'oi  6=01  s'  ;isorpm3;  ]«- '01 0=018  • 

[14]  INPUT: 'PLEASF  TYPE  ISOTOPE  CODE:  66.  60.  OR  08' 

[15]  X-ISO\J 

[16]  *(X>oISO)/ INPUT 

[17]  SO- /PC  50  [AT] 

[10]  Bl-VECBllX] 

[19]  B2-VECB2LX ] 

[20]  LAMO-VECLAMOIX] 

[21]  LAMl-VECLWllX ] 

[22]  MUO-VECMUOlX ] 

[23]  Mi/l*-V'£,C//t/l[Ar] 

[ 24  ] ' 

r//f  PARAMETERS  OF  MOLECULAR  OXYGEN  ACCORDING  TO  STEINBACH  AND  GORDY  ( 197  5 ), 
PSP.  7,  ARE  AS  FOLLOWS  FOR  ' ; ISOTEXTiX  ; ] ; ' : 

I 


[25] 

'HO  r 

' iB  0;  ' 

G5Z; 

[26] 

'd  1 = 

' ; fli ; ' 

<75  Z; 

[ 27] 

'.4  2 = 

' ;S2; ' 

cvz; 

I 

[28]  'LAM  0 = ' ;5>1M0;  ' C5Z;  ' 

[29]  'LAMl  = ' \LAMl  ; ' GHZ-, 

i 

[30]  'WO  = '-.MU 0;  ' GHZ-,  ' 

[31 ] 'MU  1 = ' ;tfyi ; ' G5Z. 

i 

[32]  TEMP- 290 

[33]  TPAfP  = ' ; TEMP',  ' REFNO  = ';  REFNO;'  ISOTOPE  = 't-TSOU] 

[34]  N INPUT-’ It2*(i27) 

[35]  -(Ar*2)/0 

[36]  sispyr-"i  + i4o 

v 


Fig.  1:  Listing  of  APL  function  PAF.AMSTEIN 
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a-E  0 d 

a boE  c-E  = °'  whel'e  J + °’  and-  . . . ( 3b) 

under  the  simplifying  assumption4-7' 11  that  B2  = 0, 

? = < J,Mj,N=j-l,s  | H | J,Mj,N=J-l,s  > 

* B°J{J-1)  + Bi  J2(  J-l) 2 + n0(j-i)  + MaJ(J-l)* 

+ 11(2/3)  - 2J/(  2J+l)^]  C Ay  + XiJ(J-lO  , 

b = < J,Mj,N=J,S  | H | J,MJ,N=J,S  > 

= B0J(J+1)  + BiJ2(  J+l) 2 - no  - Ulj(j+l) 

+ (2/3)  n Xo  + Xij(j+1)[]  , 

c * < j^Mj,N=J+1,S  I 4 j J,Mj,N=J+l,S  > 

= Bq(  J+l)  ( J+2)  + Bi(  J+l)2(  J+2)2  - uo(J+2)  - m(  J+l)(j+2)2 
+ C(2/3)  - 2(J+l)/(2J+lO  [»o*  Xi(J+l)(J+2)I]  , and 

d = < J,Mj,  N=J+1,S  | H | J,  Mj,  N^J-lj  S > 

= 2 yjj(  J+l ) (2J+1)-1  C + . ...(3c) 

In  terms  of  a,  b,  c,  and  d,  the  energy  eigenvalues  and  the  cor- 
responding eigenvector  coefficients  a,  S,  and  V in  Eq.  (i)  are 
given  by:11 


H“  - 2<B°  + 2Bl)  - 2( uo  + 2ui)  - (V3)(Xo  + 2Xr)  , 


Eo(j)  = r>. 


an  = 0xx  = yxx  = 0; 


“o  - 0,  0O  = 1,  >0  = 0; 


...(4a) 
. ..(4b) 


E+(J)  = (l/2)(a+c)  + V( l/4)(a-c)2  + d2 


‘±  " d/V(a-Et^  + a*  - S+  - 0,  T±  = a±(E+  - a)/d  . 

. . . ( 4c) 
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Eas.  { 3)  and  (4)  are  contained  within  the  APL  program  ABCD 
(lines  11-36  in  Fig.  7,  p.  24).  (This  program  computes  the  tran- 
sition matrix  elements,  as  discussed  in  a later  section. ) 

More  generally,  for  B2  / 0,  the  energy  levels  are  given  by 
Eqs.  (6)  of  Ref.  3>  as  corrected  by  Steinbach7  and  Mizushima  : 

E(  J=N=n)  = B0n(n+l)  + Bm2(n+l)2  + B2n3(n+l)3  + 2Ao/3 

+ 2Ain(n+l)/3  - ao  - fiin(n+l),  ...(5a) 

E(J=n-l)  = B0(n2-n+l)  + Bi ( n4-2n3+7n2-6n+2) 

+ B2(n6-3n5+l8n4-31n3+33n2-l8n+4)  - Ao/3 

- Ai(n2-n+4)/3  - 3..o/2  - o-i( 7n2-7n+4 )/2 

+ LB0(2n-l)  + B 1 ( 4n3-6n2+6n-2 ) 

+ B2(6n5-15n4+32n3-33n2+l8n-4)  - Ao/(2n-l) 

- ?M(7n2-7n+i+)/(6n-3)  - lXO(2n-l)/2 

- ..i(2n3-3n2+9n-4)/2H  2 

+ 4 LAo  + Xi(n2-n+l)j  2 n(  n-l)/(  2n-l)  2^J  , 

...(5b) 

E(J=n+l)  = Boi -n2+3n+3)  + Bi  ( n4+6n3+19n2+30n+l8) 

+ B2(n3+9n5+48n4+153n3+279n2+270n+108)  - Ao/3 

- \i(n2+3n+6)/3  - 3^o/2  - m( 7ns+21n+l8)/2 

- j^[>o(2n+3)  + Bi(4n3+l8n2+30n+l8) 

+ B2( 6n5+45n4+152n3+279n2+270n+108) 

- Ao/( 2n+3)  - Ai(7n2+21n+l8)/(6n+9) 

- tio(2n+3)/2  - ai(2n3+9n2+21n+l8)/2H  2 

+ 4 Q\o  + Ai  ( n2+3n+3)U  2 ( n+l) ( n+2)/( 2n+3) 2j 

...(5c) 

where  n (denoted  'n'  by  Steinbach11  and  K by  Tinkham1)  is  an 
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"effective3  (or  pseudo11)  quantum  number"  such  that  n -*  N in  the 
limit,  A/b->0.  In  the  APL  Pr°9ratn  GENGHZ,  listed  in  Fig.  2. 

Eqs.  (5)  are  used  to  compute  the  energy  levels,  ENN  * E(J-n)  -E0, 

ENNMIN1  = E( J=n-l)  = E+,  and  ENNPLU1  = E(j=n+l)  “ E_ , all  in 

gigahertz  (GHz). 

The  microwave  transition  frequencies,  f“,  are  characterized 
by  the  selection  rules, 

AK  = 0,  AJ  = t 1 (Microwave  spectrum)  (6a) 

and  are  given  by: 

f+  = E0  - E = FP  = ENN  - ENNPLU1  = ENN  - ELP  (6b) 

f“  = Eo  - E+  = FM  = ENN  - ENNMIN1  = ENN  - ELM  (6c) 

These  frequencies  are  included  in  Appendix  B ( following  the 
output  format  governed  by  the  APL  program  PRIN160O2),  together 
with  the  lower  state  energies,  ELP  and  ELM,  relative  to  the 
ground  state,  which  is  (K  = 1,  J = 0)  in  1602  and  l802  and  is 
(K  = 0,  J = l)  in  160180.  As  is  well  known, 18,21  the  homonuclear 
symmetry  of  l602  and  le02  permits  only  odd-K  states  to  exist, 
whereas  two  sets  of  states,  one  even  in  K,  the  other  odd  in  K, 

exist  in  l0Ol8O. 

The  program  GENINVCM  (Fig.  3)  converts  the  energies  to  the 
units  of  cm"1,  i.  e. , NUNN,  NUNNMIN1,  and  NUNNPLU1  are  in  cm  1 
relative  to  the  ground  state,  and  NUNPLU  and  NUNMIN  are  the  tran- 
sition frequencies  in  cm  1. 

The  submillimeter  absorption  frequencies  are  characterized 
by  the  selection  rules, 

AK  = 2,  AJ  = 0,  1 ( Submm  spectrum)  (7a) 

and  occur  as  triplets.  With  upper  states  denoted  by  a single 
prime  (K1,  J')>  and  lower  states  by  a double  prime  (KM,  JM)> 
we  label  the  transition  frequencies  as  follows: 
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[i] 
C2] 
[3] 
[4  J 

r 5] 
m 

[7] 

[8] 

[9] 

[10] 
[11  ] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[ 22] 


GENIN  V CM 

a THIS  FUNCT 10 U GENERATES  THE  ROTATIONAL  ENERGY  LEVELS 

a AND  MICROWAVE  TRANSITION  FREQUENCIES  IN  INVERSE  CENTIMETERS. 
a IT  CALLS  G ENG  HZ  AS  A SUBROUTINE,  AND  CONVERTS  TO  INVERSE  CM 
a BY  USING  AS  THE  SPEED  OF  LIGHT,  SPEEDOF LIGHT*  29.97924  58(1  . 2 ) 
a GIGAHERTZ  PER  INVERSE  CENTIMETER  ( 1973  VALUE). 
a THE  ENERGY  LEVELS  ARE  CALLED  NUNN,  NUN NM INI,  AND 

a NUNNPLU 1,  WHILE  THE  N+  AND  N~  TRANSITION  FREQUENCIES  ARE 
a CALLED  NUN  PL  U AND  NUNMIN,  RESPECTIVELY . THE  1(0)  STATE  ( FOR 
a WHICH  N~l,  J*  0)  IS  TAKEN  AS  THE  ZERO  ENERGY  LEVEL . 
a IN  THE  CASE  OF  016=018,  THE  0(1)  STATE  ( FOR  WHICH 

a N*0,  J*  1)  IS  TAKEN  AS  THE  ZERO  ENERGY  LEVEL. 

G ENG  HZ 

SPEED0FLIGHT+-  29.97924  58a  (UNITS:  GHZ  PFR  INVERSE  CM) 

+(NINPUT[1  )*Q)/ CONTINUE 
£1  0[  1 "i+ENNPLUl  [ 1 ] 

CONTIN UE : NUNN-*-(  ENN-E1  0[  1 ] ) *S?  CEDOFLIGHT 
NUNNMIN  l-(  ENNMIN 1 -El  0[1  ] ) iSPEEDO  FL IGH  T 
HUN NPL Ul+-(  EN NPLU1  -£1  0[  1 ] ) tSPEEDOFLIGHT 
NUN PLU+ NUNN- NUNNPLU1 
NUNMIN «-  NUNN -NUN  NM  IN  1 
■*(N'NPUTl  1 ]*0)/0 

HUNMINl  1 \-NUNPLUi  1 >/W////[  1 )*-NUNNMINl  [ 1 )<-NUNNPLUlt  1 J<-0 


Fig.  3*  Listing  of  APL  function  GENINVCM 
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NUF:  K«  = K“+2,  J*-1  = J..=K"  (7b) 

NUG:  K'  = K"+2,  J'  = J"  = K"+l  (?c) 

NUH:  K'  = K 1 •'  +2,  J'-l  = j.  . = K"+l  (?d) 

The  program  SUBMM02  (Fig.  4)  calculates,  in  cm"1,  the  tran- 
sition frequencies,  NUF,  NUG,  and  NUH,  as  well  as  the  respective 
lower  state  energies,  ELF,  ELG,  and  ELH  = ELG,  relative  to  the 
rotational  ground  state  (K  = 1,  j = 0 for  l602  and  leo2,  K = 0, 

J = 1 for  l60l80).  The  results  of  these  calculations  for  the 
various  sets  or  values  of  the  molecular  parameters,  B0,  Ba,  etc., 
are  included  in  Appendix  B,  together  with  the  equivalents  of 
NUF,  NUG.  and  NUH  in  GHz. 


For  the  rotational  band  lying  above  the  first  excited 
vibrational  state,  the  lower  state  energy  relative  to  the  true 
ground  state  ( for  v = 0)  is  given  by  £Gi/2  plus  the  energy  rela- 
tive to  the  rotational  ground  state,  e.  g.,  ELF.  This  considera- 
tion is  important  when  computing  the  line  strength,  and  is  taken 
into  account  in  our  calculations  (via  LINESTRENGTHROTVIB, 
APCRL5K15UPPER  and  APCRL60K1UPPER  (see  Appendix  A)). 

In  Appendix  B,  the  "laser-magnetic-resonance'  or  Raman 
lines  studied  by  Evenson,  et  al4'6  are  also  listed.  These 
lines  are  characterized  by: 

ERL:  K'  = K"+2,  J’  = K',  J*  * = K"  . ( Q) 

The  root-mean-square  deviation  between  the  Raman  lines  observed 
to  date6  and  the  various  predicted  frequencies  is  also  given  in 
Appendix  B.  These  lines  are  of  no  concern  in  calculations  of 
atmospheric  absorption  by  molecular  oxygen.17 

Line  Strengths 


The  general  expression  for  the  integrated  line  strength  of 
an  isolated  magnetic  dipole  transition  between  an  initial  state, 
(lower  state  in  absorption),  and  a final  state,  | y ^ (upper 
state  in  absorption),  taking  into  account  both  absorption  and 
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Fig.  4:  Listing  of  APL  function  SUBMM02 
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stimulated  emission  under  equilibrium  conditions  at  (Kelvin) 
temperature  T,  is  as  follows,  in  units  of  cm  1(molecule  cm-2)-1: 

S(T)  = ffe  1,0  1 < *il  ? Uf>  I 2 (e"Ei/kT  - e-Ef/kT)/Q(T),  (9a) 

where  h = Planck's  constant  = 6. 626l76( 36) • 10"27  erg  sec, 

c = vacuum  speed  of  light  = 2. 99792458( 1 . 2) • 1010  cm  sec"1, 
vo  = transition  frequency  (cm-1), 

Ei  = energy  of  initial  (lower)  state  (ergs), 

Ef  = energy  of  final  (upper)  state  (ergs), 

Q(t)  = total  partition  function  (state  sum)  (dimensionless), 
and  the  dipole  moment  operator  is  summed  over  all  degenerate 
levels  possessing  the  same  initial  and  final  state  energies.18 
Rewriting  Eq.  (9a)  in  terms  of  the  dimensionless  transition 
strengths,  l(K"  , J<  • ;K< , J • ) , 1 and  recalling  that  E - E.  = he  *0, 

S(T)  * aBohr)2  e-E"(hc/kT) 

. [J  . .-olta/Wh/ud),  (9b) 

where  gg  = free-electron  g factor  = 2*1.0011596567(35), 

^Bohr  = Bohr  "^neton  = 9. 274078( 36) * 10"21  erg  gauss"1, 
hc/k  = second  radiation  constant  = 1.438786(45)  cm  K, 

E11  “ = lower  state  energy  (cm-1), 

and  Q(T)=  ( 2K+1 ) exp(  -NUNN*  hc/kT ) + ( 2K-1)  exp(-NUNNMINl.hc/kT) 

+ ( 2K+3)  exp(  -NUNNPLU1* hc/kT )~[  . Q (T), 

QV(T)  * Cl  - exp(-AGi/2.hc/kT)C  _1  (p.123,  Ref/lS), 
since  for  each  value  of  J,  there  are  (2J  + l)  degenerate  levels.1 
(Recall  that  NUNN  = ( E0  - E )/c,  NUNNMIN1  * ( E . - E )/c 
and  NUNNPLU1  - ( E_  - Egrnd)/c,  in  the  notation  of  Eqs.  (6).11) 

The  values  adopted  for  the  fundamental  constants  are  those  of 
Ref.  13.  The  use  of  the  free-electron  g factor  is  consistent 
with  our  neglect  of  the  Zeeman  effect ; 1> *> lx* 17  an  error  in  S(t) 
of  less  than  1 part  in  104  is  incurred  thereby,  while  we  only 
require  answers  correct  to  three  decimal  places — to  be  consistent 


-16- 


T-l/306-3-14 


RIVERSIDE  RESEARCH  INSTITUTE 


with  the  accuracy  of  the  AFCRL  Atmospheric  Absorption  Line  Para- 
meters Compilation23  which  these  computations  are  intended  to 
supplement. 


In  the  microwave  region  (but  not  in  the  submillimeter  wave 
region),  « kT/hc,  so  that  one  may  replace  the  factor  in  square 
brackets  in  Eq.  (9b)  by  its  linear  approximation,  leading  to  the 
well-known  microwave  approximation: lj X1> 12 


S(T) 


± ad 

3kT 


(g, 


0 Bohr)' 2 


= e xd(  -E1  1 hc/kT) 
I(K",J"  ?K',J')  vo  Q (T) 

(Microwave  approximation)  ...(9c) 


(Note  that  in  the  optical  region,  the  factor  in  square  brackets 
in  Eq.  (9b)  is  approximately  unity,  i.  e.,  the  induced  emission 
term  may  be  neglected  at  300k  and  below,  as  in  Eq.  (3)  °f  Ref.  23.) 

Choosing  T = T = 296K,  the  standard  temperature  of  the 
s 

AFCRL  compilation, 23  the  thermal  energy  expressed  in  wavenumbers 
is  kT/hc  = kTs/hc  = ( 296 )/( 1.438786(45) ) = 205.729(6)  cm"1,  so 
that  Eq.  (9b)  becomes: 

S(Ts=296K)  * 1.4353-10"22  I(K",J"  ;K»,  J')  *0  exp(-E"/205.729) 

• Cl  - exp(  - >'o/205.729)II  /Q(Ts=296k)  . 

...(9d) 


Since  we  desire  our  computed  line  strengths  to  be  consistent 
with  those  of  AFCRL, 23  we  must  multiply  S(Tg)  by  the  relative 
isotopic  abundance  of  the  molecule  (our  Table  I,  or  Table  3 °f 
Ref.  23),  which  is  0.99519  for  the  dominant  isotope,  le02.  De- 
noting the  resulting  quantity  by  Sm(Tg),  whose  units  are  cm 
per  molecule  of  mixed  oxygen  per  cm2,  we  obtain  an  expression 
like  Eq.  (9d)  with  the  numerical  factor  in  front  replaced  by 
the  quantity,  1.4283967* 10"22  (retaining  extra  digits  to  avoid 
round-off  error)  which  appears  in  line  22  of  LI NE STR ENGTH 2 
(Fig.  5). 

The  rotational  state  sum  q(t)/Q  (T)  is  evaluated  by  summing 
Eq.  (10)  through  K = 79  in  the  program  ISOSTATESUM  (Fig.  6),  and 
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Fig.  5:  Listing  of  APL  function  LI NESTRENGTH2 
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V ISOSTATESUM-tKT\N  INPUT \TERMliTERM2-tTERM3\TERM 

" THIS  PROGRAM  COMPUTES  THE  ROTATIONAL  PARTITION  FUNCTION 

[2]  n OR  STATE  SUM  FOR  MOLECULAR  OXYGEN.  IT  IS  BASED  ON  THE  ENERGY 

L3|  a LEVELS  OF  STEINBACH  AND  GORDY.  THE  PROGRAM  SETS  TEMP  EQUAL 

[4]  " TO  KELVIN  TEMPERATURE  TYPED  IN  RESPONSE  TO  'QUAD'  (J)  PROMPT. 

AT  296 K,  THE  LASTJTERM  INCLUDED  IN  THE  STATE  SUM  FOR  016=016 
LGJ  A CONTRIBUTES  ONLY  8.0  65T17;  FOR  016=018,  THE  LAST  TERM  IS  8.40£“l6. 
[7]  N INPUT*-  ((2*l40)-l) 

[31  +{X*2)/ SKIP 

[9]  NINPUT*-{  18  0 ) -1 

[10]  SKIP : GEN INVCM 

[11]  'PLEASE  TYPE  IN  KELVIN  TEMPERATURE : ' 

[12]  TEMP*-.] 

[13]  KT*-  2 0 5 . 7 29  * ( TEMP  *296) 

[14]  TERM  1 «-  ( 1 + 2 *//  INPUT  )*(*((  -NUNN  )i  KT  ) ) 

[15]  TERM2*-(~\  +2*N  INPUT  ) *(  *(  ( -NUNNMIN 1 )i'KT)  ) 

Tl  G]  TERM3*-{  3 +2*N INPUT  )*(*((  -NUNNPLU 1 ) iKT)  ) 

[17]  ->{X*2)/ADD 

[131  TERM  1 [ 1 ~\*-TERM2ll  1*-0-,7ERM3[_l  ]«- ( 7;i7?73  [ 1 ] ) *3 
[13]  ADD-.T  ERM*-  TERM 1 +TERM  2 * TER  M3 

[20]  SUM*-* /TERM 

[21]  CLASSICAL*-^  1.  S* KT*SPEEDOFLIGHT)  *50 

[22]  'ROTATIONAL  STA  TF.  SUM  = ' ; SUM ; 1 ( EXACT  VALUE)' 

[23]  *(X-2 )/BYPASS 

[24]  'CLASSICAL  APPROXIMATION  = ( 3i  2 )*(  KT  i B 0)  *S  PEEDO  FLIC  NT  = ' -.CLASSICAL 
[2  5]  CONTINUE  :RA  TIO*-SUM  iCLASSICAI, 

[26]  'RATIO  = RATIO 

[27]  'TEMPERATURE  = '\TEMP\'K\  ISOTOPE  = • ; ISOTEXTlX ; ] ; ' . ' 

[ 20  ] -*•  0 

[ 29  ] BYPASS -.CLASS  I CAL*- 2*  CLASSICAL 

[30]  'CLASSICAL  APPROXIMATION  = ( 3 *K7 iB 0)*SPEEDOFLIGHT  = ' -.CLASSICAL 

[31]  -*•(  A'  = 2 )/ CONTINUE 
7 


Fig.  6:  Listing  of  APL  function  ISOSTATESUM 
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compared  with  the  classical  values, x*  l2' 18  3kT/2hcBo  for  1602 
and  l802,  and  3^T/hcBo  for  l60l80.  (The  state  sum  contains 
twice  as  many  rotational  states  in  a heteronuclear  diatomic  mole- 
cule as  in  a homonuclear  one.18)  Appendix  D lists  some  results. 

The  vibrational  partition  function  QV(T)  has  been  approxi- 
mated by  unity  at  T = T = 296K.  This  is  consistent  with  the 
AFCRL  approximation  of  Q (T  ) * 1.000  (Table  2,  Ref.  23).  More 
precisely,  in  the  harmonic  oscillator  approximation  given  by 
Eq.  (10),  CQV(TS=296K)H  -1  = 1 - exp(-AGi/2-hc/kTs) 

= 1 - exp, -1356. 378/203. 7?9) 

* 1 - 5.1818-10"4  * 0.999^82,  (11) 

leading  to  an  overestimate  of  5.2  parts  in  104  in  the  values  for 
Sm(Tg)  as  given  in  this  report.  Thus,  the  third  decimal  place 
in  the  line  strengths  tabulated  in  this  document  will  be  in 
error  by  one  unit,  at  most. 

We  now  turn  to  the  calculation  of  the  normalized  transition 
matrix  elements  (squared),  I(K' », J» • ;K», J' ),  i.  e.,  of  the  quan- 
tities F,  G,  H,  PP,  and  MM  in  lines  27-31  of  LINESTRENGTH2 . 

Transition  Matrix  Elements 

Steinbach11  has  shown  that  the  normalized  magnetic -dipole 
transition  matrix  elements,  both  those  for  the  microwave  "fine 
structure"  (or  spin-reorientation18)  transitions  and  those  for 
the  "forbidden"  submillimeter  rotational  transitions,  may  be 
written  in  terms  of  the  Wigner  6-j  symbols20' 24  and  the  eigen- 
vector coefficients,  a and  ?,  appearing  in  Eqs.  (l)  and  (4). 

With  n'(J')  denoting  the  final  state,  n(j)  the  initial  state, 
i.  e.,  employing  n instead  of  K and  dropping  the  double  primes, 
the  results  are  displayed  in  Table  II.  These  expressions  may 
be  simplified  further:  First,  we  note  that  for  fixed  J ^ 0, 
o_  = and  a+  = -y_  (c£.  Tables  8,  9>  and  10  of  Ref.  11).  To 
prove  this  formally,  let  us  first  note  that  d > 0 from  Eq.  (3c). 
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Table  II:  Normalized  Magnetic  Dipole  Transition  Matrix  Elements 


Transition 
n(  J)  - n'  ( J'  ) 


a 

Normalized  Matrix  Element 
I Of  l"l*i>  I 2/(gseflBohr)2 


Pine  Structure  Transitions: 
n(  n-1)  n(  n) : 

n(  n+l)  -•  n(  n)  : 

Rotational  Transitions: 
n(n)  **  n+2  (n+l)  : 

n(n+l)  -*  n+2  (n+l):  6(2n+3)2 

n(n+l)  -•  n+2  (n+2)  : 


eian-lXSn+Dr2!^  J J} 

6(2n+l)  ( 2n+3)  { n+f  J l} 

6(  2n+l)  ( 2n+3)  n£2  { * "f1  "}2 

v*{^,r?}+ VfUil11!8 

6(2n+3)(2n+5)\2{nl1 


a 


After  Table  2,  Ref.  11.  The  quantities 
Wigner’s  6-j  symbols  (see  Ref.  24). 


ji  J2  j 3 \ 
Jl  J 2 J3  ( 


are 


T-l/306-3-14 


-21- 


RIVERSIDE  RESEARCH  INSTITUTE 


Letting  (c  - a)/2d  = £ in  Eq.  (4c),  we  find  that 


x±  s (E±  - a)/d  = i + yj  1 + 

so  that  x+  x_  = £2  - 1 - = -i,  and  x_  = -l/x+.  Letting 

x+  = tan  0,  Eq.  (4c)  then  states  that  a+  = cos  0,  7+  = a+x+  = sin  0. 

Since  x_  = -l/x+  = -ctn  0,  Eq.  (4c)  also  statefi  that  a = sin  0, 

= a_x_  = -sin  0 ctn  0 = -cos  0.  Thus,  a = y+J  a+  = -y  f qed. * 


Next,  using  Eq.  (C.36)  of  Ref.  24  to  evaluate  the  6-j 
symbols,  we  obtain: 


n 


< 1 n n l 2 = ( n+l ) . ( 1 n n I 2 _ „ 

' 11-1  1 1 I 6n(  2n+l)  ' I n+l  111  ~ 6(  n+l) (2n+l)  ' 

) 1 n+l  n I 2 _ n J 1 n+lnl  / (n+2)  ~~ 

' n 11’  6(  n+l)  ( 2n+l)  ' I n+l  1 1 J y6(n+i){  2n+3)  ; 

j 1 n+l  n+2 1 = / (n+l)  J 1 n+2  n+2  I 2 = (n+3) 

’ n+l  1 1 f \6(  n+2)  ( 2n+3)  ' I n+l  1 if  6(n+2)(2n+5) 


Inserting  these  results  into  the  expressions  in  Table  II,  noting 


that 


n+2 


L "n+l 


n+l 
I n+2  J 


- ( 2n+3) 2( n+l)/( n+2) , and  replacing 
n by  K once  again,  we  obtain  I(K’  ’,J’  ’ ; K ’ , J ’ ) = j_5_  I H I ^ >1  2 


T 


0 


Bohr 


y 


MM  = I ( K,  K-l  ;K,  K)  = C >+(K-l)^  (2K-1)(K+1)/K  ~ (i2) 

PP  = l(K,K+l;K,K)  = C o_(K+l)^]2  (2K+3)K/(K+l)  (13) 

F = I ( K,  K ; K+2,  K+l ) = H 7_(  K+l  )I]2  ( 2K+3)  • K/(  K+l ) (l4) 

G = I ( K, K+l ;K+2, K+l ) = C *_( K+l )H2( 2K+3) • £ ajK+l)^ 

• ( 2K+3) 2( K+l)/(K+2)  (15) 

H = I ( K,  K+l  ;K+2,  K+2 ) = £ M K+l  )J2  ( 2K+3)  • ( K+3) /( K+2)  (l6) 

where  the  J-dependence  of  the  alphas  and  gammas  in  Eqs.  (12) -(16) 

is  indicated  explicitly  by  their  arguments  in  parentheses. 


From  Eq.  (l),  it  is  intuitively  clear  that  0 represents  the 
small  but  non-zero)  angle  in  Hilbert  space  through  which  | ib.y 
are  rotated  away  from  their  Hund’s  case  (b)  approximations.  - 
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The  equivalents  of  these  six  equations  are  contained  in  lines 
4l,  42,  45,  and  46  of  the  program  ABCD  (Fig.  7).  (To  avoid 
possible  division  by  zero,  the  dyadic  APL  function  X REPLACE  Y 
is  used.  It  is  listed  following  ABCD  in  Fig.  7.  It  replaces 
zero  values  of  J and  K by  a small  number  during  arithmetic  opera- 
tions (ABCD  line  10)  and,  at  the  end  of  the  calculations  (line 
4l),  it  undoes  the  effects  of  this  replacement.) 

It  is  clear  from  Eqs.  (l4)-(l6)  that  the  submillimeter 
transitions  would  be  strictly  forbidden  were  actually  zero 
(pure  Hund's  case  (b)  coupling)  instead  of  merely  small.11  On 
the  other  hand,  (and  a_)  are  nearly  equal  to  their  case-(b) 
values  of  unity;  thus,  the  microwave  transitions  (Eqs.  (12)-(13))^ 
are  negligibly  affected  by  the  deviation  from  case  (b)  coupling. 
The  various  types  of  allowed  transitions1*  7’  8* 11  are  shown  in 
the  energy  level  diagrams  of  Fig.  8. 

The  existence  of  a common  factor  in  Eqs.  (l4)— (l6)  not  only 
simplifies  the  calculation  of  the  strengths  of  the  "forbidden” 
lines,  but  also  enables  us  to  exhibit  the  asymptotic  behavior  of 

their  relative  strengths,  since  a_(K+l)  -•  1 as  K -°°,  e*> 

F:G:H~1:4:1.  This  is  in  agreement  with  the  correspon- 
dence-principle argument  of  Gebbie,  et  al.8  On  the  other  hand, 
in  their  final  results  (Table  VII,  column  6,  and  Eq.  (63a)), 
Tinkham  and  Strandberg1  (T-S)  incorrectly  and  " unphysical ly" 
predicted  ”F"-type  transitions  to  be  much  weaker  than  H -type 
transitions,8  an  error  carried  over  into  the  Russian  literature. 
Nevertheless,  the  intermediate  results  of  T-S  (Eqs.  (60)  and 
Table  V)  were  essentially  correct,1  as  pointed  out  by  Gebbie, 
et  al.8  Indeed,  one  easily  deduces  that  F ; H = K(K+2)/(K+l)  (K+3) 
both  from  our  Eqs.  (I4),(l6)  as  well  as  from  the  fifth  and 
eighth  of  Eqs.  (60)  of  T-S.1  Thus,  except  for  I(K,KjK+2,K+l), 
our  results  for  the  transition  matrix  elements,  given  in  Appen- 
dix B,  are  in  close  agreement  with  Table  V of  T-S. 

To  within  the  tabulated  accuracy  and  limited  spectral 


T- 1/306-3-1^ 


-23- 


RIVERSIDE  RESEARCH  INSTITUTE 


CO  to 

XC  Si 
Jf  St 
a:  to 
xc  3: 
ti  to 
>o 
Cl  to 

fis  to 

x:  ^ 

• CN 

• ft:  H 4* 

tO  Si  V 

*<  X<  • ••  • 

^ atj  o 
to  tt  • ^ 
ai  %to  V S: 

o <r.  ^ a. 

•sq  s o 

to  Si  St  " xt 

xt  s,  *?  a. 

a:  to  tt  • 
a,  % Si  ••  3-  • 
so  ..  to  z-.  ^ Si 
•orn;  Vca 
Si  ci  • ft.  a. 
to 

Ci  to  xr  ..  «t  c> 

o t*:  «h  o 

n Si  a i . 

ft:  «*  :<  a.  to 

xt  to  a;  • ^ xt 

^ Si  V * a: 
d a-  w a.  to 


at 

c 

a: 

S; 

o 

tj 

• • 

O 

5: 

o 

•• 

tj 

X 

Si 

to 

'll 


a c 

t Qo 

to  »i 

- . -,  a.  ci  ci 

to  xt  •-1  ft.  ♦ Cl  to 

a-  v t» 


an r o to  v 

Ci  ta  ~ 

• Si  Si 
xa  Cl  II 
a: 


. tN 

ta  o,  oj  v ^ c£ 
o ci  ••  a, 

tt  ’o  to  » to 

a,  % to  - ♦ oa  to  % 

si  *o  ^ K s.  li, 
'-O  ft:  O •••Stto 
t a - a<:  a.  Si  Si 

ctr  t«i  • v-  to 

Si  a.  ft:  ><sn  c- 

a.  ci  w ii  ti*  h 

^ ti,  s ti  x a; 


■o> 

X 

vH 

5: 

■f 

o 

5: 

•-3 


>:  ~ 
w ^ 

♦ 2r 

^ + 

o 

r-1  1!^) 

1 5; 


- • 

r.;  • V- 

.. 

1 

Si  to  ft- 

*s 

CM 

* 

z—x 

ti)  |~1  xc 

Vo 

4 • Co 

o 

z<^ 

3:  to. 

* 3E  n 

Or 

n to  to 

••  xt 

5: 

X 

ft;  ft:  Si 

xi  * as 

* 

Xt  Si  Cl 

t« 

♦ a-  ‘-i 

4- 

3: 

•• 

a.  a. 

o 

^ xt  ft) 

Z-X 

C3l 

to  Si 

• ta  s. 

t) 

z«*v 

^3 

•• 

to  - Cl 

to 

CM 

4- 

5:0 

C3 

x— z • 

^3 

4 

O 

*• 

Si  *>  . 

> 

xx  a.  ^ 

cx* 

Z*“S 

» 

a* 

"t  to 

" xt  ft; 

U: 

4: 

•• 

to  to  ft: 

a:  ?;  -t 

ft: 

H 

ix3 

Im 

a;  a. 

1 

•• 

S,  *1 

&; 

to  xj; 

o 

X 

r- 

■*  :1 

O 

a:  ta  to 

* 

z*^ 

o 

•-5  Si  • 

ft: 

xt  >-4 

i 

X 

cn 

o 

d to  ft; 

co 

• x4 

U: 

•*• 

ft; 

to  U: 

- O xr 

H 

XwZ 

CM 

•• 

►4  • 'a 

xt  to 

»«z 

XwZ 

"o 

it  ft;  t*i 

• 

*-*  t-x 

ot 

X 

XwZ 

> 

to 

Z-K 

♦■  to  3; 

Si  .» 

4- 

• o 

• 

V 03  xj 

to> 

cn 

•• 
f X 

• ft; 

a. 

X— ^ 

CM 

c* 

z— v 

^ Cl 

CM  ^ 

4 

^ • 

K 

o 

O 


>: 


o 

* 


fti  a.  ^ 

^ ^3  ^ 

z«^ 

a:  >•  Si 

II  ^ 03 

CM 

Or 

Si  p 

4 

O 

to  t.. 

5:  • ^ 

z-S 

X 

Q O 

Ox  Si 

tr~i 

rl 

*5  Cl 

••  Cl  «t 

1 

> 

Xt  Si  =5 

z->  s,  S. 

*"3 

Cl 

V Si  xt 

vz 

• 5;  *~i 

• *n 

w 

4* 

CUtOK 

V O i 

X 

O 

a:  xc 

.•  Cl  ft; 

•o, 

5: 

• Si  Si 

•X  "t  xt  to 

x ■ — ■ 

'J  Cl 

♦•  Sr. 

XI  Ct 

M5os 

Ci  x 

o 

>: 

ft: 

Co 

as 

Co 

as 

Co 

o 

ft.  ^ 

Co  ^ 

Uj  Z*“S 

P ^ 

Co  *-» 
* 

ft: 

o 

e-,  * 

Cj  CM 
Co  4- 

^ o 
as  *; 

r*j  n; 

O *0 


Co 


x 

CM 

4 ^ Q n 
•h  ^ >*.  •*• 
w ^ ^ CM 


a:  to 

s,  tj 
=~  x: 

+ so 

v a, 

..  Co 

*<  ft:  zi 

• *-i 

z—  O ♦■ 
n J1! 
•*■  I w 
t<  to  x 
-rt'l 

* * A. 

‘-s  * 


* 

^ 

^ X CM 
H rt  ♦ 
I 03  ^ 
w w 
w 4-  X 

X z—x  z—x 

^ O'  fH 
X x 4- 
O O ^ 

03  03  w 

>*✓  >*✓  I 

♦ ± 2 

ni  XX'y 


4- 

•H  Qr 
4 w 
X 

w *H 

* 5: 

CN  ^ 
w iO 
I w 
^ 4- 

m o 

4 J* 
CM  ^ 
w kO 
w 'w' 
w X 
4-  ^ 

Z*S 

O'1-) 
cy  x 

X CM 

^ 4- 

cy  th 

X w 
r-i  •*• 
03  z-s 

w m 
4-  • 

z-x  O 

ty  * 
cy  Of 

* w 
O X 
CQ  CM 


tO|  il 


to  ^ 

Co  4 

ro  -« 

0: 

Co  5: 

^ * 
t*3  CN 
Is  «xZ 

r.3  4 

'o  o 

Co  5: 

»wZ 

c«)  x 

rt:  cn 

t-i  — 

t 

to  -. 

Si 

CS  n 
a.  ai 
3:  * 

8C; 

o 
i to 

0 ^ 

% X 

CN 

to  — 

1 * 

Ci  an 
^ + 
a;  ^ 
vj  o 

C O to 


>0 

o 

a: 

to 

to 

Ci 

•J 

Xt 

tx 


<t 

xi 

tl 

to 

a. 

to 

g 

o 

•o 

Sig 

Cl 

0 to 

Cr.  S; 

1 Si 

z“*x 

CM  » 

in  •♦•  Co 

H 

o R.  ^ ♦ 

* n Si  ti 

-i  w to  n 

cn  + 5; 

* •a.  ta  a 

ciia  ti  n 
<•  fa  ci  tj 
••  Q •• 

tjt  •*  o ta  * 

1 ^ o ta 

<N  Ck;  Co  n 

> * •*•  * Si 
tj  tj  —>  Si  to 

a a:  <n  oa 

HI  "■  "5  •• 
^ ’O  S O 
UwB,  i 
♦ >-z  xr  ta 

^ » vz  ix 

> Si  + 
tj  c>  > 
a<  o ti, 
r a 5i 


§ 

to 

N 

■5 

o 

% 


5 

05 

>i  O 
xl  Q 
>.  — ■ 
O * 

Q 

• 

*5.  + 

O 

" ? 


in  m 

^ XC 

ta 
ii  .• 


a: 

a. 

vj 

■ <c 


a 

fe  z-* 

V 
a:  a 
Si  ix 
•»  * 
i o 
xc  o 
'-0  Si  Q 
to  V •”* 

Si  05  V 

hO  v 

to  a * 

«”x  a: 

tk;  V xt 
S Q Si 
5.xil 
* 03 
O * 

3 

0.  Si 
xa  '-*5 

C ^ 03 


o o 

* * 

z*x  z*\ 
CM  CM 

* • 

C^C3I 


o • 

65  ^ 

Co  Ux 
O 


Si 

• 

• xt 

o<  a: 

s3 

s. 

S3 

n 

Si 

to 

Si  xt 

Cl 

Si 

a.  to 

»t 

a 3: 

*~l  X( 

Si 

a:  to 

Ci 

ft. 

O to 

s. 

S,  a: 

Ci 

Si 

Cl 

to 

to  to 

to 

a-  si 

a> 

to 

C3  % 

is. 

ax  o 

Xt 

£>  ^ 

3: 

a.  to 

to  a« 

a: 

to  to 

to 

ft:  3: 

Si 

ft:  si 

a: 

o c> 

4 

to 



to 

ft:  a: 

to 

to  Si 

fti 

a. 

^5 

Ci  to 

D5 

ft:  te 

xc 

ft.  Cl 

Si 

to 

ft: 

ax 

xt 

Sx  to 

£x 

5x- 

Si  XC 

xft 

S3  3: 

xt 

>1  3- 

S:  ta 

ax  n 

xt  Q 

Si  to  a: 

3:  s5 

o a:  a. 

3:  ta 

ax  Si  xo 

XT 

_ 5:  o « 

xc  *■  ft: 

ft;  ax  '— 

Si  Z-.  xt 

Xt  Xt  V 

Ink 

CM  CM 

« 4 

z*^«  S 

i i 


c^lC>| 

ii 

X X* 

a:  a: 
a.  a. 
«o  *o 


Q 03  Q tO 

to  -0  xx  +■  _ «• 
ft;  to  * —'  * 3: 
o Cl  3:  ^ s-  3: 

har  Q ^ ^ 

Co  ^ ft;  w 3;  ca 

ftl  5x  CU  * 5 

^ o 3:  xt  q 

- xt  - ta  a: 

to  - ••  Si  ••  xt 

Co  ^ C*3  z-s 

to  Si  a.  03  a,  . 
?Jt  J ••  x(  to 
Si  Si  a zi  a ^ 
tai  cn  to  a,  - - ^ 
* * » o ^ x 
^Si  o n g ta  03 


xciitta  , 
• i sa 
-»  ^ o o 
aian; 
to  Ci  ta 

w w to  • 
X 

3:  a: 

^ Si  to 
ta  ft; 
tk.  Si  Si 
•o  xt 


C3 

a: 

to 

a. 

a, 

Xt 


os  ta  • 
• xt  to 
ta  3:  3 

3 % ? 

fi  xt  a. 

'•3  ta  oj 

05  i «( 

w 4- 
4-  z-s  bg 

XI  V ft: 
V ,3  f*< 
a - 

T >i 

to  5 

05  ta  c 


' .n.i.  .n.nnnnnnnnnr ii— inr- 1 1— inni— i 

^7  . 77  7OrtlNnj‘n't"v®l7'OHrvn*1'Mi,r't“it,Ort<xnjini0N 

TitNm^ioiCtioDCCTiTirtrti-iiHiHWxxi-itNfMfMCNtNCNCNfMfMtNriririctinrinm 

*— * SJ  UJ  U LJ  l_J  UJ  LJ  txJ  k i I I I i I i i a i a ■ , a , a a . , . . . 


Fig.  7:  Listings  of  the  APL  functions  ABCD  and  X REPLACE  Y 
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coverage  of  Steinbach's  thesis, 11  all  of  our  transition  matrix 
elements  in  Appendix  B agree  with  those  of  his  Tables  11-13. 
Steinbach  has  also  noted11  the  discrepancy  of  these  results 
with  those  tabulated  by  Tinkham  and  Strandberg. 1 We  have  been 
unable  to  ascertain  whether  or  not  our  results  also  agree  with 
those  of  Gebbie,  et  al,  who  only  give  integrated  line  strengths 
in  units  of  10  6 cm  2 atm  1.8  In  Appendix  E,  we  list  our  results 
for  integrated  line  strengths  in  these  units,  taken  to  refer  to 
2?3K,  so  that  l(cm-atm)STp  = 2 . 68675M 84) • 1019  molecules/cm2. 23 
The  results  in  Appendix  E exhibit  discrepancies  of  as  much  as 
a factor  of  two  with  the  results  of  Gebbie,  et  al.,8  but  inas- 
much as  no  temperature  is  specified  in  Ref.  8,  the  precise  mag- 
nitude of  the  discrepancies  is  difficult  to  ascertain. 

Line  Widths 

For  molecular  oxygen  at  pressures  characteristic  of  alti- 
tudes below  40  km,  pressure-broadening  predominates. 26  Near 
an  isolated  transition  frequency,  pressure-broadened  molecular 
absorption  is  well-approximated  by  a Lorentz  line  shape  with 
half-width  • = *c*p  at  half-maximum  points,  where  p is  the  pres- 
sure. 23,27  (Bold-face  alphas  are  used  to  distinguish  them  from 
the  quantities  appearing  in  Eqs.  (l),  (4),  (13),  and  Table  II.) 
Then  the  absorption  at  the  peak  of  the  spectral  line  will  be 
given  by  S0  = s(t)/jt«,  where  S(t)  is  given  by  Eq.  (9)  ,Z3’ Z7*  and 
will  fall  off  away  from  the  peak  in  accordance  with  whatever 
lineshape  is  most  suited  to  the  particular  spectral  region:28  Van 
Vleck-Weisskopf  in  the  microwave  region, 12  "kinetic"  in  the 
submillimeter  region, 29,30  and  Lorentz  or  "super-Lorentz"  in  the 
infrared  and  visible  regions.31 

Of  the  line  parameters  required  to  perform  atmospheric  trans- 
mission calculations, 23  the  linewidths  for  oxygen  are  the  least 
accurately  characterized  ones,  unfortunately.  The  experimental 
values  of  «o  (in  MHz/torr)  measured  for  l60160  in  the  microwave 

region  by  several  investigators  between  1952  and  1968  (at  300K) 

* S(T)  actually  should  be  replaced  by  Sm(T).  (See  page  17. ) 
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have  been  summarized  conveniently  in  a review  article  by  P.  H. 
Krupenie.32  These  experimental  values,  converted  to  cm  1/atm 
and  rounded  to  two  significant  figures,  nave  been  employed  in 
the  line  parameter  compilations  for  molecular  oxygen  contained 
in  this  report,  (l  MHz/torr  = ( 760 ) /( 29979* 2458)  - 0.02535  cm  1 
per  atm.)  To  the  degree  of  uncertainty  (as  much  as  25$)  inher- 
ent in  the  data, 26,33  no  correction  has  been  made  for  a possible 
isotope  effect  in  the  line  widths;  instead,  a crudely  linear 
interpolation  on  K was  made  from  Krupenie' s Table  34,32  to  des- 
cribe the  widths  expected  for  the  even-K  states  in  160;L60  (see 
Fig.  8),  and  the  tabulated  values  themselves32  were  used  for  the 
odd-K  states  in  all  three  isotopes  considered. 

Other  than  the  relatively  crude  results  in  Ref.  8,  no  data 
are  available  in  the  literature  concerning  the  widths  of  the 
submillimeter  lines  of  molecular  oxygen.  If  one  assumes  that 
the  line  widths  depend  on  the  value  of  K in  the  lower  state, 
and  one  notes  that  "G"  and  "H"  transitions  with  K'  ' = n possess 
the  same  lower  states  as  the  microwave  "+"  transitions  with  the 
same  n (see  Fig.  8),  then  one  may  safely  assume  that  these  two 
types  of  submillimeter  transition  possess  the  same  line  widths 
as  the  associated  microwave  transition.  Somewhat  arbitrari- 

ly, the  line  widths  for  "F"  transitions  have  been  taken  as  the 
arithmetic  mean  of  the  linewidths  of  the  and  transitions 

with  the  same  value  of  n. 

An  alternative  procedure  of  even  less  accuracy  would  have 
been  to  use  the  electron-optical  data34-36  to  deduce  the  sub- 
millimeter line  widths.  However,  unlike  the  K-dependence  found 
for  the  microwave  line  widths,  the  J-dependence  exhibited  by  the 
line  widths  in  the  visible  region34-3 6 is  unlikely  to  be  correct 
for  the  submillimeter  lines. 

For  lines  corresponding  to  higher  values  of  K than  have 
been  measured  in  the  microwave  region,  a lower  limit  of  0.032 
cm-1/atm  has  been  assumed.  The  line  widths  are  entered  in 
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lines  27-30  of  AFCRLFILE  (Fig.  9),  lines  29-31  of  AFCRL60K15 
(Fig.  10),  and  similarly  in  the  other  "AFCRL"-format-setting 
programs  listed  (by  name  and  short  comment  only)  in  Appendix  A. 

The  determination  of  line  widths  near  room  temperature  and 
1 atm  for  molecular  oxygen  is  complicated  by  the  demonstrated 
breakdown  of  the  independent-line  approximation  in  the  microwave 
region. 2°’ 33  The  microwave  line  widths  have  been  treated  in 
the  past  by  assuming  «o  to  be  pressure-dependent,  3 7} 38  and  more 
recently  by  means  of  several  interacting-line  theories.39  42 
Since  the  microwave-line-width  situation  in  molecular  oxygen 
is  under  current  active  study  at  several  institutions, 43  it  was 
decided  to  use  the  Krupenie  compilation32  to  obtain  interim 
(±10  to  2556)  values  for  the  line  widths,  pending  further  results 
expected  during  the  next  year.* 

Line  Parameters  in  "AFCRL  Format" 

McClatchey,  et  al23  have  described  their  compilation  of 
molecular  spectroscopic  parameters  for  seven  infrared-active 
molecules  which  occur  naturally  in  the  terrestrial  atmosphere. 

As  of  1 April  1975, 44  all  of  the  submillimeter  lines  of  oxygen 
were  absent  from  the  AFCRL  computer  tape,  and  the  microwave 
lines  were  all  represented  as  having  the  line  width,  0.060  cm  1 
per  atm.  It  was  in  order  to  supply  the  data  on  molecular  oxygen 
which  was  missing  (and  is  required  as  input  to  the  RRI  SLAM 
program45)  that  the  calculations  described  in  this  Research 
Note  were  undertaken.  Having  performed  the  calculations  des- 
cribed in  the  preceding  sections,  we  therefore  made  use^  of  whe 
APL  output-formatting  program  AFMT  (described  in  Ref.  46-)  to 
cast  our  results  into  a form  compatible  with  the  AFCRL  format.23 
This  was  accomplished  by  means  of  the  eight  APL  programs  whose 
names  begin  with  "AFCRL"  in  the  defined-function  list  in  Appendix 
A.  Two  of  these  (for  l60160  in  the  vibrational  ground  state)  are 

* Since  S(T)  has  been  calculated  fairly  precisely  in  this  report, 
measurements  of  peak  absorption  in  the  submillimeter  region  should 
yield  line-width  values  which  can  help  clarify  the  microwave  line- 
width  situation  (because  of  the  relationships  discussed  above). 


T- 1/306-3-14 


-29- 


RIVERSIDE  RESEARCH  INSTITUTE 


<0 

<0 


tn 

r- 

n 

X 


X 

£ 

£ 


*0 


a: 

o 


C0 

co 


£ 

a 

% 

to 

aE 

tt;  to 
tj  to 

fiS 

2^ 

h • ^ 

•*  X Ofe 

X O =a 
X % 

.•  U] 
o ^ 

! Or 
« to 
i a: 
a. 


• to  • 
mo'' 
H to  to  «B 

ti 

to  to  % 

Cl.  X t_> 

V X 

a a a 

oc  to 
to  a.  a. 

X 

H to  to 

►sat*. 

• o<J  <o 
to  to  >- 
H a a 

>c  a o • 
t«  a % * 

a X * 

» a x 
a:  to 

to  to  a.  • 
is; 

Otil  * 

a G *, 

a o a • 
X to  - 

'""EL 

« to  N 
a:  * a. 
to  a;  % • 


to 

X 

2 

* 

<o 

<0 


to 

to 

a 

a 

to 

a. 

to 


“•e 


te8. 


a 


JS  to 
X X 
..  *, 
a 

*»  £ 

■ • ri 

x to 


a a x 


•s  * 

to  •< 
••  a: 
to  a 
••  o 
a a; 
a a. 

x to 

••  a 

to  a 

a h 

a 

a. 

a 

a 

a 

a 


a 
a 
r< 
Q 

__  to 

a a 

oh,  a 
a 

a a a • 

X a a - 
a a ^ 
a 

* a a • 

a a h * 
a a G 
a ••  , 
a to  a a 
a a a a 
*■•  o a a 
a * . a 
v a 
a a h a 
or  a m h 
a a 

• a • 
a . % 

Ha  a 

yt'jt 

*"  o a a 
a a a % 
o m a 
a a a . 
a ^ a to 
a a a 
o •Ha 
a a a a 
H * •* 
a v a *t 
a a a h 


a a 
a 
ats 
an 
q a 
a q 
a 

a n 

• r~ 

a a 

2 a a a 

xaag 


a a - _ _ 

a a ••  h ° ° 

a a - ••  a ° a 

_.  o » n.  o . 

a ■ a~a 

a a H a a a 

to  sax 

lfet.ni  n 
a a a % -<  «**  •-' 
% ^ a a * *•»  * 

a:  a a * a 

« « - - a aa 


a 

a 

a "-i 
± x 

'Ou 

<n  a 

• a 

**a 

X ♦ 

* ^ 
tn  *4 

w-  I 
U-l^ 

^ X 
^ x 
••  tn 
ft.  w 
hOu 
♦ 

X X *5 


x x 

nnn  i_j  uj  m 

h x x a a x 

uuun  a a a< 

a a a x a a a 

a a a i—i  ♦ ♦ a 

fe  t ( ferinfe 

i m m a tN  ..  4 

n w n * * i rn 

n*  i i 

W <-»  ->  X X X 
* X X X * * * 
(•)  n x x o n n 

w rt  M <•>  W w *-* 

■ , w w w l— < t— I > — l 

^uu>-<  ^3  ^ ^ 
X *0  <0  <0  t«3  ^ ^3 


04 

tn 

t4  3 O 
X tn  • 
o o o 


X 

£ 


O O CM 

n 

tn 

• 

•H  * o 
a tn  • 

12 

o o o 

» 

H 

O O 

X 

tn 

• 

CM  10  o 
3 tn  • 

CM 

X 

o o o 

• 

• • 

O O in 
tn 

MOO 

c 

3 tn  • 

o 

o o o 

FI 

• • 

• 

O O 

10 

tn 

X 

CM  CO  O 

• 

3 m • 

o o o 

12 

• • 

m 

O O »0 

H 

tn 

X 

m co  o 
3 tn  • 

• 

CM 

o o o 

X 

• • 

• 

. O O N0 

a 

fM 

r~ 

<r 


a 

a 


ot  a n 

aQ  » o'  o 

a ^ rl  • 

a a ° ° ° 
* a • 
p n ° 
a a 
a » 
a a * 

a d o 

a 

a a-  o 

£g- 
5*2 

si; 

to  ^ 

C*  vt.  * 
to  o 
<0 


! *.  • 
; o o 

l m 


fXfeCO 

05; 

oh 

OCi 
Co  ^ 

4*  VI 

x x^ 

X ♦ « 


CNN  * 

o 
C0  ° 

Qo  00 
S 

to  o 
to  • 
a° 

C ^ 


O k0 

tn 
O'  o 
tn  • 
o o 
• 

o <0 
tn 
*4  o 

3 • ^ 

O o CM 

• tn 
0(00 

tn  • 
a-  o O 
o*o 
. o ^ 
O 

10  b« 
a tn  a 
o o x 

• • tn 
o o w 

♦ 

H 0 <0 

& tn  tn 
o o i 


X 

*H 


h3 

X 


*4 


V 


U tx 

X X 

..  H 

x a- 

UJ 


n 

tn 


CM  CM 

♦ ♦ ^ 
X X 


4 

♦ 4*  w 

XX* 


o o w 
• • • 

iii 


4 

m X X 

X X • 
— . W X 

3>  a a 

*•»  a *» 

h ^ ^ ^ 
s, 

luuu 

a a a a 
h •»  x ►, 
a > ♦ ♦ 
♦ a a a 
x ►>  x •» 


a u> 

- a 
• a 

m ^ • 

nN  tjj  r> 

X-  • 
. . o 
X-  -i 

• a a 
x-  • 

4 — m 
a a • 
X nn  O 

w a •• 
x a a 

I— I * » 

iiS 

Stih 


X o 

i — . — 

a a ■ 


Sj 

a ..i-i 
WLJH 

a a H 
* a •• 
a •<  >-< 

< o ci 
w a 
Si  aa 

x — ■ a 

I«4  — ' — ' 

a * a a 

4 i-i a 
a *-i  * a 
fe  (11-10 
fe  XH  I 
ouitxia 
..  d ••  •• 
M feuh 

a •<  a a 

— i a a a h 

♦ <j  fto 


X 

X 

Cto 

X 

X 

o 

Cto 


X 

o 

<0 

2 

r«- 


CO 

to 

X 

s 

t-5 

O 

co 


8 

£ 


a 

o 

to 
a 
I tar 
. a 
, a 

is 


a a 

x a 
r<  a X 
f VI  O 
X X 
* *•» 

X ♦ « 


nnnn  — nnn— 1/1 
„p— .nnnnnnn  On«nSVlOr.»0'  H(in  SNMVN(vnnnO  n nn  n n " *,5.5.^  ^ 

Htins  v>or.«<r  x h h r*  ••  cv  incv  « ^ ^ ^ uj  lj  uu  vj 


Fig.  9:  Listing  of  APL  function  AFCRLFILE 


-30- 


t-1/306-3-11^ 


RIVERSIDE  RESEARCH  INSTI1UTE 


FJ 


VO 

VO 


• t/J  • 

G ft;  Or 
N ft;  G 

* * 
tV  * G 

>1  5 


• ft;  tv  I II 

ft:  t ^ ^ id 

O * ft;  O 
SMijIijo 
* >-3  * kj  to 

to  S a 4-  ft: 

% 

•-1 
•-3 


ft: 

t. 

Cl 

ft;  sc 
•t  tc 

••  £-i 

E-.  ^ 

>< 

ft;  Ui 
E-.  G 
••  N 
H »-3 
E-. 


‘ ft; 


O 

o 

VO 

N 

G 

r. 


<N 

o 

o 

VO 

E-. 

ft: 

H 

ft; 

ft. 


ft;  G 


••-JO* 

ft] 

ft: 

h]  Uj  w 

G 

to 

to 

t.  ^3  uj 

Jr 

'w’ 

* £t  O • 
t.  O ^ • 

VO 

to 

N 

to  * 

10 

ft: 

S-. 

^ ft:  * 

w 

» 

ft: 

to  ft:  ft,  .to 

VO 

:*  ft  -to 

H 

Oht-  N 

O 

-3  O -»  O 

II 

ft.  ft: 

ci  0 ft:  . n. 

H 

* to  - to 

O M 

Tt-ft- 

ft: 

* ft:  ^ O 

ft;  to 

itift  *— 

P P 

ft, 

*•4 

G 

g 

ftj 


v. 

ttj 

G 

ft; 

G 

ft, 

ft: 

o 

G 

G 

ft; 


*H  G 


G 

ft: 

n: 

to 


to  • to 
ft*  ft: 
> G G G 
i H N 
i G ^ 
I to  ft: 


<7> 

CO 

O 

O 

CD 

CO 

o 

o 


a 

o 


O' 

CO 

O OD 
• CO 

o o 


ft:  rl 


CM 

CO 

o 

• 

o 

in 

CO 

o 

• 

o 

in 

co 

O 


<N 

N 


CM 

S 


CN 


FJ 


o 

C 


to 


CM 


>< 


»-3 


SO 

ft: 


..  • 

G O 

h - 

O 

0 

H ^ * 

N 

O 

• 

4—1 

a * 

ft,  ft. 

* 

t- 

■ - 

* * ta  tv 

• 

r- 

rj 

G 

* 0 

53 

• 

O 

•j  •j 

* * tt;  to  O 

O 

m 

*•»  ^ 

..  * 

ft: 

ft  • 

to 

sc  * 

t.  SC  to  sc 

4-3 

0 

••  H 

•3  ft: 

ft,  ft: 

G - 

N 

0 H 

to  * ••* 

ft: 

H 

* 

• 

^ ft; 

*»  ^ 

O E-. 

r,  S 

* 

tv  1 to  to 

G 

ft. 

X 

0 

0 

N ft. 

• ■ Of 

ft; 

ft. 

to  * 

O -~l  tc, 

H 

ft: 

0 

1—1 

O to 

to  ft: 

to  • 

O 

H sc 

to  * to  V 

i-3 

^4 

0 

• 

in 

G ^ 

* ft 

ft;  H 

ft;  - 

to 

W ta 

ft. 

5: 

G 

0 

CO 

Ui~ 

••ft. 

ft; 

ftj  ^ 

to 

S:  • 

0 * ^ * tv 

G ^ 

0 

n 

• • •* 

=3  , 

ftj  • 

5* 

f-. 

•c  * 

to  sc  to  0 

^5  fM 

N ft; 

ft3 

CM 

• 

• 

ro  0 

s ft 

ft:  ft> 

G • 

H 

* * 

*1  SC  O tt;  to 

f«;  ft: 

»-3  G 

N 

J- 

O 

4-N  fO 

M • 

..  cv 

53 

^ to  • 

ta  to 

0 to 

ft:  p 

ft,  H 

ft 

O 

• 

4—1  a • 

* 

ft:  53 

G 

*~i  — 

O * 

SC  + IV  to  to 

ft,  H 

G *-J 

ft; 

• 

CD 

^ O 

G r 

• • 

to 

•• 

O sc 

* * 

to  % to  t. 

"-0  »-3 

P ft- 

ft. 

O 

CO 

t4 

• • r » 

G to 

ft:  ft; 

ft;  ft; 

to 

* * 

Si  VI  > to  C: 

•*•  ft. 

ft;  0 

G 

0 

• ft. 

ro  ro  q 

••  ti 

53  i-3 

ft:  to 

0 0 

Qf  to  0 tv  So 

^ O 

ft;  eft 

ft; 

r- 

• 

CM 

m 

M ^ 

4-3  to 

^ P 

* ^ 

tO  .H 

to  ■•  *-S  ti  3c 

r-s  Q 

ft,  h3 

x 

a 

0 

CO 

w ro 

l—l  CN  <0 

G H 

ft:  ft: 

ft; 

0 * 

N O 

iv  i-i  <c  * 

m ft; 

G ft3 

0 

O 

^■S 

• • 

G 

/• 

to  ^ 

tj  0 

* -J 

to  „ to  Kj 

4—i  ft; 

ft. 

ft: 

• 

in 

• 

^■s  m 

• •l—l  *H 

to 

53  G 

^ ft: 

* 

tv 

1—1  tv  S>  to 

ft,  ft. 

^ G 

0 

CO 

O 

• 

tH  ft. 

ro  p CN 

Of  ft: 

N N 

to  II 

* 

SirtVtl 

ft,  G 

r-4  ^ 

0 

a 

X 

4-  - G 

M ft;  •• 

G 

• • sc 

* * . OJ  * 

w •*■ 

N ^4 

G 

x 

• 

ft:  1 • 

1— 1 ^ 4—1 

ft:  ft: 

G • 

• r~ 

to  to 

h V to 

4 ^ 

4—1  rn 

ft: 

X 

0 

4 

w • CO 

G G G 

••  to 

G • 

* 

* ~ 

n.  1 

S:  * 0 tv 

f“l 

a.  n 

G 

0 

G 

w • * 

■n 

G w ft; 

M O 

E*  G 

* . 

ft;  0 Jft 

to  * ^ C~i 

N 

•3  1-1 

• 

as  ft. 

a-  0 

M 

W Q ^ 

m a: 

* “• 

• T 

60  * 5)  * 

* Sc  c * 

1 t— 1 4—1 

r- 1 to  3t 

G 

0 

CO 

a 

^ 4*  r-l 

CO 

w G 

••  n. 

**•  O 

ft  ^ 

^ »-3 

ft:  O ft; 

^ to  S;  CO  i; 

^ ft:  ft. 

N ^ S3 

G 

0 

X 

to 

i—l 

ft:  W W 

M 

to  E-. 

P * 

5)  Sto  * 

1 lO  * 

O 

O 

M ft,  G 

l_l  + to 

O G 

in 

• 

CM 

4—1 

* W . 

£-4 

ft,  w a 

G 

••  to 

O 

N O 

ft  O Nft, 

iitovs; 

a 

O 

a 

X w 4 

O G 

0 

0 

w» 

G 

a au1 

G 

< ro  w 

G 

in  h 

ft;  to 

^ . 

* to  * 

tv 

* a ta  *: 

^■s 

a 

(N  1 n t/]  h i 

4-3  £-1 

• 

♦ 

w — t 

O 

rl  4 

rj  ftj 

ft;  H 

*r  fti 

p 

ta  •• 

•3  Cl  SC 

ft. 

ft. 

wr-i  ^ 

4 • ^ 

^ G 

0 

*• 

O' 

4 

• a.  0 

G 

E-t  CN  ro 

G 

»-3  to 

0 

G 

* to 

Gft,ftft;c3ft«ft,ft, 

1 1 ^ l 

^ N 

CO 

CM 

t-4  G 

tN  ft,  r,  33 

4 

••  n 

4 

0 

G ,* 

N ft; 

ft  ft:  so  ^ 

* * -c  tv  * 

a tv 

a 

53  M ^ 

^ M M 

ft,  ^ 

in 

0 1 

— ■ a*  * 

G 4— l CM 

g 

lO 

G ft; 

n.  * 

te  ft: 

* * to  * 

K 

a 

X 

^ X M 

M X X 

ft,  G 

X 

• 

tv 

4 

w a - •* 

•• 

tiG  «• 

*-3 

Eo 

S:  ^ * Q 

rl 

^ H 

to  tv  to 

CM 

X 

CM 

"MX 

X CM  CN 

iH 

a G 

0 

0 

W’ 

* 

G 

»«*  * 

CN 

4 G 

h 

ft; 

ft,  ^ 

0*0* 

II 

• ^ 

tv  « to  a. 

V— * 

CM 

w* 

ft,  w CN 

CM  w w 

4- 

Vi  H 

• 

• 

• 

N 

5< 

•wV  h 

ft. 

G G 

g 

0 

O G 

O H i:  ■« 

K 

ft: 

to  tv  to  sc  tv 

w 

4 

»H  Q 4— 1 W 

wu  u 

M M O 

0 

1 3 

S*s 

4 

4 4 * * 

4 

O 

G G ft; 

N 

N— ' 

0 P 

53 

4 

4 

C3  53  »— » 

*-»  S3  *-3 

4 

s— |» 

4 

i 

4 

4 

4-3 

•3  * to  0 N Q 

^ wt 

G 

^ c 

c c 

c c 

c c 

t 

t G 

1 v 1 j i 

tnft;Mk4ft:^;^;ft;ft3 

M 

4 c 

4 

* 

3 

’■»  to  fo  to  N ^3 

G 4 G 

G 

> 

ro 

m 

m 

rn 

m 

nn  n 

nn  n 

ri 

m m 

ro 

m 

r-i 

ro 

*"o  nn  ro 

ro 

ro 

ro  ro  ro 

ro 

ro 

r*o  ro 

m ro 

ro  ro 

m 

ro  0 

h n cn  ^ in 

vD 

r- 

00 

<T» 

O H M 

cn  ft-  m 

m 

00 

O' 

O 

fH 

CM  CO 

^ m in 

00 

O' 

O rl  CN 

n 

ri 

tN  CO 

X 

vD  1*. 

00 

CT'  rH 

H H H H H 

H 

rH 

H 

CM  CM  CM 

CM  CM  CM 

CM 

CM  CM 

CM 

co 

CO 

CO 

CO 

CO  co  ro  CO 

CO 

CO 

* * * 

X 

i—J 

l— 14—1 

L— J l_J 

V— i l— i 

i—J 

4—1  4— J 

4—1  4—1  4—1  4—1  4—1 

4—1 

4—1 

4— 1 

4— J 

4— i 4—1  4— i 

4— i 4— J 4— J 

l—l 

l—l  4—1 

4—1 

4—1 

4—1 

4—1 

4—1 

4—1  4— J L— J l—l 

4—1 

l—l 

4—1  4—1  4—1 

4— J 

Fig,  10:  Listing  of  APL  function  AFCRL60K15 
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iisted  in  Figs.  9 and  10.  The  format  specifications  in  line  37 

o RLFILE  (Fig.  9)  and  line  35  of  AFCRL6OKI5  (Fig.  10)  are 
equivalent  to  the  FORTRAN  format: 

strength  WIDTH  E'»  v ' T • v, 

( 10.5,  Eio.3,  P5.3,  F10.3,3X,I1,3X,I2,1X,I2,6X, 

II,’  3X,  K,  IX,  12,'  IX,  A3, 5x,  Diaf  2X,  lt%x,  S) 

This  format  is  in  accordance  with  the  one  given  on  page  6 of 

Places  (nTT  thT  T frS,UenCy  13  9lven  to  «ve  decimal 
Places  (F10.5)  instead  of  three  (F10.3). 

In  Appendix  F,  the  contents  of  the  eiaht 

by  APCRLP1LE,  APCRLK2  ArcRLfinKin  generated 

j AFCRL60K35  are  listed.  (There  arp 

two  fries  each  for  the  three  isotopes  (‘VM,  ‘V.0  V*ol 

in  the  vibrational  ground  state  (v  - 0)  and  in’the  flrst 

excited  vibrational  state  (v  = c h flrst 

lines  th^  * 1 ' f°r  thG  submillimeter 

lines,  the  other  for  the  microwave  lines.)  As  indicated  by  the 

heading  at  the  beginning  of  each  file  listing,  the  line  para- 
meters were  generated  by  using  Steinbach  and  Gordy-  for  the 

tional  PTetCrS  (B°'  Bl'  B*'  1*0.  Ui)  in  the  vibra- 

i onal  ground  state,  and  Albritton,  et  ,i«  for  those  of 

® V £tate*  The  function  LINESTRENGTH2  (Fig  5)  Dr 
s abbreviated  form,  LINESTRENGTH22  (see  Appendix  A),  ms  used 

compute  the  line  strengths  at  296k.  (for  convenience  the 
vibrational  partition  function  was  ignored,  i.  set  ^ to 

unity,  as  discussed  above.)  The  line  widths  were  derived  from 

i:;p:r:„:ompiiatW2  in  the — “ - - 

coi  P°r  clarityi  at  the  toP  of  each  page  in  Appendix  P,  the 
columns  have  been  labelled  with  the  appropriate  headings: 

SSZ-il  (cm-Amoiecule  cm-));  line 

St  (cm  atm  );  lower  state  energy  EH  (cm’1),  upper  state 

quantum  numbers.  V'  jt  i HP  te 

v„  a . „ " ^ and  lower  state  quantum  numbers, 

' — ' n — ; sh^thand  line  identification  (K+  or  K-  for 
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the  microwave  lines;  SF,  SG,  or  SH  for  the  submillimeter  lines); 
month  and  last  digit  of  year  in  date  of  computation  (75  = July 
1975.)*  isotope  code  (66  = l60160,  68  = 160l80,  88  = 1801B0) ; and 
molecular  constituent  (j  = oxygen). 

The  eight  files  listed  in  Appendix  F were  merged  into  a 
single  file,  ordered  by  frequency  by  means  of  a FORTRAN  program 
(SRWRTMERGE)  which  employed  the  XDS  Sigma  9 MERGE  processor. 

A threshold  criterion  was  then  applied  to  the  resulting  merged 
file  to  select  only  the  most  "significant"  lines:  only  lines 
whose  strengths  exceeded  3-7  * 10  30  cm  1 per  molecule  cm  2 
(the  "Existing  Intensity  Minimum  at  T = 296K"  given  for  02  in 
Table  3 of  McClatchey,  et  al23)  were  to  be  retained.  (The  zero 
frequency  "lines"  of  zero  strength  in  the  listings  for  160l80 
in  Appendix  F — artifacts  arising  from  the  manner  in  which  the 
APL  programs  were  executed — are  removed  at  this  stage. ) The 
resulting  file,  OXYGENEXIST,  is  listed  in  Appendix  G.  (This 
file  was  created  by  means  of  another  FORTRAN  program,  NOSIG, 
which  allows  one  to  remove  "cards"  whose  strengths  lie  below 
a user-specified  level  from  any  card-image  file  written  in 
"AFCRL  format".) 

The  file  OXYGENEXIST  in  Appendix  G contains  the  parameters  of 
all  "fine-structure"  and  "rotational"  lines  of  the  molecules, 
l60160,  l60l80,  and  l80l80,  which  can  be  expected  to  be  of 
significance  in  atmospheric  absorption  problems.  A deck  of 
punched  cards  containing  the  card-image  records  in  this  file 
was  sent  to  AFCRL  on  28  Tuly  1975* 

Summary  and  Comments 

Calculations  performed  at  RRI  have  been  described  which 
give  current  best  estimates  of  the  line  parameters  of  all  sig- 
nificant microwave  and  submillimeter  wave  absorption  lines  of 
the  molecular  oxygen  isotopes,  l60l60,  l60l80,  and  l80l80,  at 
zero  magnetic  field.  The  transition  frequency  (given  to  the 
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nearest  .00001  cm  1 - 0.3  MHz)  and  energy  of  the  lower  state 
(quoted  to  the  nearest  .001  cm  1)  have  absolute  accuracies  on 
the  order  of  .001  cm  1 as  judged  oy  the  degree  of  consistency 
between  measured ’~7  and  predicted  values  (Appendix  B) . These 
results  are  certainly  in  agreement  with  recent  low-resolution 
stratospheric  Fourier-transform  spectroscopic  measurements4 6 
of  the  submillimeter  emission  by  the  earth's  atmosphere  between 
35  and  200  cm  1 . The  rms  deviation  between  our  calculated 
frequencies  and  the  twenty-seven  observed  lines  of  l60i60  short 
of  200  cm  1 is  .05  cm  1 , which  is  within  the  .07  cm  1 (unapo- 
dized)  resolution  of  the  measurements. 

Integrated  line  strengths  at  296k,  accurate  to  within  one 
unit  in  the  third  decimal  place,  have  been  calculated;  the  new 
results  are  in  apparent  disagreement  (Appendix  E)  with  earlier 
published  line  strengths8  for  i60l60,  by  a factor  of  approximate- 
ly two.  We  note  that  these  earlier  line  strengths, 8 together 
with  the  strength  data  on  water  and  ozone  lines  from  the  AFCRL 
compilation, 23  have  been  used  recently  to  deduce  the  temperature 
of  the  cosmic  background  radiation, 47  as  well  as  to  determine 
the  water  content  of  the  stratosphere. 48' 49  A closer  look  at 
the  data  analyses  may  be  warranted  in  both  cases46-49  to  deter- 
mine whether  the  deduced  results  are  sign.*  ficantly  affected  by 
the  change  in  oxygen  line  strengths  to  the  new  values  we  report. 

The  least  certain  parameters  are  the  line  widths,  for  which 
we  chose  interim  values,  based  on  Ref.  32.  The  line  widths  in 
the  microwave  region  are  currently  under  intensive  study41-43 
so  that  we  should  be  able  to  improve  on  the  line  width  estimates 
during  the  next  year.  We  pointed  out  that  direct  study  of  the 
submillimeter  line  widths  could  help  to  clarify  some  of  the 
issues  attending  the  more  complicated  problem  of  the  microwave 
line  widths.  The  line  widths  which  we  give  here  must  be  viewed 
as  uncertain  to  +10  to  25^. 
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APPENDIX  A 


Brief  Description  of  APL  Notation  and  Output  Formats 

The  intent  of  this  Appendix  is  to  enable  the  reader  unfamil- 
iar with  APL  to  verify  that  the  algebraic  expressions  given  in 
this  report  have  been  implemented  properly  within  the  APL  func- 
tions” (programs) , and  to  explain  the  APL  output  formats  in 
Appendices  B and  E. 

APL50”52  is  an  interactive,  time-sharing-oriented,  program- 
ing language  which  is  especially  suited  to  working  with  arrays 
in  one  or  more  dimensions  with  a minimum  of  looping  or  branching. 
It  uses  a special  symbol  set,51'52  many  of  whose  elements  have 
intuitively  clear  meanings  (e.  g.,  the  parenthesis  pairs,  ( and 
),  or  the  binary  operators,  +,-,♦,  etc.),  but  the  same  symbol 
may  play  different  roles,  depending  on  the  context  (e.  g.,  0 or 
p).  The  following  rules  are  of  principal  importance  in  following 
the  APL  programs  listed  in  this  report: 

1.  Operations  are  performed  from  right  to  left,  with  no 
hierarchy  between  addition,  multiplication,  etc.  In  most  of 
the  programs  listed,  parentheses  have  been  used  liberally  to 
make  the  resulting  expressions  look  more  like  ordinary  algebra. 
Examples:  ♦ means  "multii-1Y  by  reciprocal  of  ev^rythina  to  the 
right";  -/  means  "sum  the  alternating  series  resulting  from  the 
placement  of  a binary  minus  between  all  elements  of  the  array 
sequence  beyond  the  slash  (/)'  • 

2.  Operators  can  be  unary  or  binary  ( *B  means  "take  recip- 
rocal of  B" ; A+B  means  "A  divided  by  B"  where  A and  B are  arrays 
of  the  same  "rank"  (dimensions)  or  scalar);  functions  can  be 
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mladic  (e.  g.,  GENGHZ),  monadic,  or  dyadic  (e.  g.,  REPLACE,  or 
the  intrinsic'- 1,52  function  £FMT)  . 

3.  Substitution  is  indicated  by  an  arrow  pointing  left; 

M - 5 (equivalent  to  FORTRAN:  M = 5)  is  not  the  same  as  M = 5 
in  APL,  the  latter  (equivalent  to  FORTRAN:  M.EQ.5)  having  0 or 
1 V no  or  yes  ) as  a result.  (if  M is  an  array,  the  logical 
relations  return  arrays  of  the  same  dimensions,  filled  with 
zeros  and  ones.) 

4.  APL  has  two  kinds  of  minus  signs.  Negative  numbers, 

such  as  in  the  exponents  in  Appendix  C or  in  front  of  the  ele- 
ments of  VECB1  in  PARAMSTEIN,  are  preceded  by  the  unary,  upper- 
case minus,  ; negation  or  subtraction  is  indicated  by  the 
binary,  ordinary-looking,  minus,  -.  Example:  If  A is  "5  4 0 "3, 

the  result  of  the  operation  -A  is:  5 ~4  0 3. 

5.  While  standard  texts50-52  should  be  referred  to  for 
the  meanings  and  syntax  of  the  special  APL  symbols,  the  follow- 
ing occur  sufficiently  often  in  our  programs  to  merit  discussion: 

, ( comma J used  as  a binary  operator  means  "concatenate": 
012,  01  produces  the  result  01201.  (See  lines  27-30  of 
AFCRLFILE, ) Used  as  a unary  operator  (as  in  line  34  of  A'ORL60- 
K15),  it  means  "ravel",  or  string  out  as  a vector.  (The  special 
symbol  \ means  "transpose".) 

* ( semi— colon)  is  used  for  separation  of  independent  ex- 
pressions, as  well  as  (in  headers — line  o)  to  indicate  local 
variables  defined  only  within  the  program  (as  opposed  to  global 
variables  such  as  SPEEDOFLIGHT,  defined  "forever"). 

* (asterisk)  denotes  exponentiation:  A*B  means  AB  (binary 
operator);  *B  means  eB  (unary  operator). 

1 (vertical)  denotes  absolute  value:  | A means  I A |. 

' (x°ta)  is  the  "index"  operator,  with  unity  as  the  origin: 

1 xi5  produces  the  result  “l*(l  2 3 4 5)  or  "l  -2  “3  “4  “5 
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p (rho)  used  monadically  is  the  dimension  operator:  p A re- 
turns an  integer  vector  with  the  dimensions  of  A.  Example: 
p(2  ~1  5)  returns  3*  (The  integer  vector  has  one  element  in 
this  case  since  A was  a vector.)  P (rho)  used  dyadically  is 
the  reshape  operator:  Ap  B returns  an  array  whose  dimensions 
are  given  by  the  left  (vector)  argument  A,  and  whose  (initial) 
values  are  given  by  the  right  argument  B.  If  B does  not  contain 
sufficient  values,  it  is  recycled.  Example:  2 5 p 3 returns  the 
2-by-5  array,  33333. 

3 3 3 3 3 

t and  + are  the  "take'  and  "drop"  operators.  Examples:  Let 
A be  the  array  332403.  Then  3*A  yields  332;  3+A  yields 
4 0 3.  See  also  lines  34-39  of  the  function  ABCD  (Fig.  7). 

[ and  ] (brackets)  enclose  line  numbers  in  program  listings 
and  also  indicate  array  indices.  Examples:  If  A +-  2 3 1 0, 

then  A[3]  is  1;  if  M is  a two-dimensional  matrix,  M[5;]  picks 
out  the  fifth  row  (all  columns)  of  M. 

6.  Branching  is  indicated  by  the  arrow  pointing  right,  and 
the  syntax  is  illustrated  by  the  example:  -*•(  x^2)/NEXT  which  is 
equivalent  to  the  FORTRAN  statement  IF  (X.NE.2.0)  GO  TO  NEXT. 

7.  The  APL  system  carries  out  all  calculations  to  approxi- 
mately 16  significant  digits.  The  results  are  displayed  rounded 
of£  to  the  number  of  digits  specified  by  a system  command, 
)DIGITS.  Even  after  output  display  is  called  for,  the  full  16- 
digit-accuracy  is  retained  in  the  system  memory. 

8.  The  "default"  output  format  for  the  vector  quantities 
in  Appendix  B is  to  be  understood  as  a list,  in  one-to-one  cor- 
respondence with  the  array  NINPUT  given  at  the  start  of  each 
page.  In  Appendix  E,  the  "take  5"  operation  has  been  applied 
before  printing  out  results  relating  to  Ref.  8,  so  that  these 
arrays  refer  to  the  first  five  values  of  NINPUT.  (Similar  com- 
ments hold  for  the  experimental  values  of  ERL  in  Appendix  B. ) 
Note  that  trailing  zeros  are  suppressed  on  display. 
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”BLE  A'1:  0F  APL  FUNCTIONS  CR EAT ED_FOR_TH?S "report 


function  NAMF 


COMMENT 


PAGE 


ABCD 

ABCOLP 

ABCDUPPER 

afcrlfile 

AFCRLK2 

AFCRLK3 

AFCRL5K1 SUPPER 
AFCRL6UK1UPPER 
AFCRL60K1 5 
AFCRLG0K25 
AFCRL6  UK3  5 

COMPARAMAN 

FREIIL  1ST 

GBBSTRtNGTH 


G ENG HZ 
GENINVCM 


ISOSTATESUM 

ISOSTATESUMLP 

LINESTP.ENGTHROTVI 
L INESTRENGTH2 
L INESTRENGTH22 

PARAMETERS 

PARAMETER SLP 
PAR AM STEIN 


PA  RAM  STE  I NLP 
PARAMSTEINU PPER 

PRINTGHZ 
PR  INTI  NVCM 
PR INTSTATESUM 
PR INTSTRENGTH2 
PR  INT SUBMM02 
PR  INT  GOO 2 


RFPLACE 

SUBMM02 


ucfc.I?  calculate  transition  matrix  elements 

uPDcmN  0F  ABCD  USE3  ,T°R  L,NE  PRIkTrER  OUTPUT 
VERSION  OF  ABCD  FOR  V = 1 STATES  OF  016=016 


FOR 

FOR 


016=016 
016=01 8 
01  8=01 8 
016=016 


FOR  V = 0 
(V  = 0) 

CV  = 0) 
FOR  V = 1 


LISTS  LINE  PARAMETERS  OF 
VERSION  OF  AFCRLFILE  FOR 
VERSION  OF  AFCRLFILF 
VERSION  OF  AFCRLFILE 

GHZ'  L,NES  0F  016=016  FOR  V = 1 
60  GMZ'  L,NES  0F  016=016  FOR  V = 0 
vpdc !2N  °F  AFCPl 6 OKI  5 FOR  016=018  (V  = 0) 
VERSION  OF  AFCRL60K15  FOR  018=018  (V  = 0) 

COMPUTES  RMS  DIFFERENCES  BETWEEN  PREDICTED 
AND  OBSERVED  LASER  'RAMAN'  LINES 

S!f!IL^UTfJUT  °F  PARA”STEINLP,  SUBMM02 , 

PR  I NT  6 00  2,  AND  COMPARAMAN  (SEE  APPENDIX  B) 

a*“°upares 

WITH  RESULTS  OF  REF.  8 (SEE  APPFNDIX  E) 

™!!5c«TES  Ef,ERGY  LEVELS  IN  GIGAHERTZ 
CONVERTS  GENGHZ  OUTPUT  TO  INV  CM  A M3  GFNFratfc 
MICROWAVE  TRANSITION  FREOUENCIES  (ALSO  INV  CM) 

VFPP  ,0NAL  PART  IT  ION  FUNCTION 

VERSION  OF  ISOSTATESUM  FOR  LINF  PRINTER  OUTPUT 

B VERSION  OF  L INESTR ENGTH2  FOR  016=016.  V = 1 
COMPUTES  LINE  STRENGTHS  FOR  ALL  OXYGEN  LINES 
COMPUTES  LINE  STRENGTHS  OF  SUBMILLIMETER  L^NES 

^ANY^VrU  LVAM0'  LAM1'  MU0'  MU1  ACCORDING 
v;2oc^L°NE  0F  PEFS-  1 THROUGH  6 FOR  016=016 

GWF  °Bf)°Fw?AR!l9FTERS  F°R  LINE  PRI‘JTER  OUTPUT 

LAMO,  LAM1,  MUO,  MU  1 ACCORDING 
7 F°R  ANY  OF  THREE  ISOTOPES;  SETS 
NINPUT  AND  TEMPERATURE 

ilpRcmM  °F  PARAMSTEIN  FOR  LINE  PRINTER  OUTPUT 
VERSION  OF  PARAMETERS  FOR  016=016,  V = 1 (REF.  5) 

PRINTS  OUTPUT  OF  GENGHZ 

PRINTS  OUTPUT  OF  GENINVCM 

PRINTS  STATE  SUM  VS.  TEMPERATURE 
PRINTS  OUTPUT  OF  L INESTR ENG TH 2 

PRINTS  OUTPUT  OF  SUBMM02 

E NERg'ie^6 S a c cd F TRANS  IT  10 NS,  LOWER  STATE 
ENERGIES,  AND  LASER  'RAMAN'  LINES  (IN  GHZ) 

AND  PRINTS  RESULTS  (SEE  APPENDIX  B)  ' 

BINARY  FUNCTION  WHICH  REPLACES  ’ZERO'  VALUES 

^C.U!;ATES  su»hillimete:r  transition  frequencies 

AND  LOWER  STATE  ENERGIES  (IN  INVERSE  CM) 


(see  append ix  d; 


LISTINGS 

(MARKED 


24-25 

N/A 

N/A 

30 
N/A 
N/A 
N/A 
N/A 

31 
N/A 
N/A 


N/A 

N/A 

N/A 

12 

13 

19 

N/A 

N/A 

18 

N/A 


73 

N/A 


8 

N/A 

74 

N/A 

N/A 

N/A 

N/A 

M/A 


N/A 

25 

15 


OF  APL 
N/A  IN 


FUNCTIONS  MENTIONED  HERE  BUT  NOT  APPEARING  EXPLICITLY 
THIS  COLUMN)  ARE  AVAILABLE  FROM  THE  AUTHOR  SSfl  REO^ST. 
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APPENDIX  B 


Molecular  "Rotational"  Parameters.  Millimeter  and  Submillimeter 
Wave  Transition  Matrices.  Transition  Frequencies,  and  Lower 
State  Energies  in  the  Ground  (32^~.  y = q)  State  of  l60l60. 

— - — — lB0le0.  and  in  the  First  Excited  Vibrational 
S£.3te  LY  = l) — of  1 60<2 i According  to  References  1 through  7 


For  each  of  the  ten  sets  of  parameters  studied,  there  are 
three  pages  of  computer  printout.  The  first  of  the  three  pages 
gives  the  rotational  parameters,  BO,  Bl,  B2,  LAMO,  LAM1,  MUG, 
and  MU1  in  GHz  and  in  cm  1 , the  Reference  number,  isotope  code, 
the  array  of  K-values,  and  the  corresponding  arrays  of  normalized 
transition  matrices.  The  second  page  repeats  the  Reference  num- 
ber and  array  of  K-values  (NINPUT),  followed  by  the  submillimeter 
frequencies  (cm  and  GHz)  and  lower  state  energies  (cm-1)  with 
respect  to  the  ground  state.  The  third  page  repeats  the  Refer- 
ence number  and  NINPUT,  followed  by  the  microwave  transition 
frequencies  and  lower  state  energies  with  respect  to  the  ground 
state  (GHz)  and  the  predicted  and  observed  "laser-magnetic-reso- 
nance" lines6  and  their  rms  difference  in  MHz  and  cm-1. 

The  sequence  of  the  ten  sets  of  parameters  is  as  follows: 
Ground  vibrational  state  of  l60i60  (v'  = v' ' = 0),  in  order  of 
Reference  number  (1  through  7);  Ground  states  of  l60i80,  l8ol8o 
^Ref.  7);  first  excited  vibrational  state  (v'  = v' ' =1)  of 
160160  (Ref.  5).  (In  the  last  case,  AG i/2  is  not  included  in 
the  lower  state  energy,  as  discussed  in  the  text  on  p.  14.) 
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APPENDIX  C 


ADDITIONAL  APL  PARAMETER -SETTING  PROGRAMS 


The  parameter-setting  program  PARAMSTEIN  (Fig.  1,  p.  8) 
gave  B0,  Bx,  B2,  \o,  Ax,  ,l0,  and  ux  for  le0le0f  l6ol80,  and 
O O,  according  to  Steinbach  and  Gordy, 7 for  the  vibrational 
ground  state,  v = 0.  These  seven  parameters  (for  l60l60  only) 
according  to  Refs.  1 through  6 (for  v = 0)  are  contained  in  the 
APL  program  PARAMETERS  (Fig,  C-l) . The  APL  program  containing 
the  parameters  for  the  upper  vibrational  state:,  v = 1,  in  le02, 
according  to  Ref.  5,  is  PARAMSTEI NUPPER  (Fig.  C-2) . In  this 
program,  LGi/s  is  represented  as  GV. 
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V PARAMETERS ; X ; VECBO-.VECBl  iVECB 2;  /PCZ.4.V0;  VECLAMl ; VECMU0-.VECMU1 : REFT  EXT  \ -.REFTEXTi  -.BEE TEXT}  ; REFTEX  T*  ; 
RE‘-,7EXT6-.REFTEX7 X 

' ' 1 • 'HEFNO ' /S  4 V INTEGER  FROM  1 3*0  6.  CORRESPONDING  TO  ONE 

[?!  • JP  .’#£  S«  REFERENCES  LISTED  IN  ’REF?  EXT • . 

[’1  /.VS0*43. 1 929  43.100589  43.100518  4 3 . 1 00  51  8,  ( 1 . 4377  0»*SPEEDO FLIC HT ) ,4  3.  1 004608 

! 4 1 />TS1**0. 0001471499  *0.00014  *0.0001449628  *0 . 000144  96  , ( *4  . (4  E‘ ~6”SPE  EDO  FLIC  NT ) , *0 . 0001 4 5? 

[5]  /PCS’*  0 0 *1  . 57  i*”l  0 *1.7  S*1  0 0 0 

[’J  / ECLAM  0*59.  501  57  59.501346  59.501342  59 . 501 34  2,  ( 1 . 5«1 . 3239  •SPEEDOFL 1GNT) . 59 . 501  3 42 

t’J  VSCLX’ll-i.  678S*5  5.  84  5F*  5 5.847P*5  5.847f*5  0 5 847S*5 

.'£T.V7  0**0.2  52  67  *0.2525917  *0.2525865  *0 . 2 52  58  6 5.  ( *0 . 0084  36*  SPEEDOF  LIGHT  ) . *0 . 2 52  58  6 5 
[9]  /.-’CV  71*0  *2.45  5£~7  *2.4640*7  *2.464f*7  0 *2.464E*7 

110]  REFTEXTI-  '? INKHAM  AND  STRANDBERG' 

[11]  REFTEXTI-  'WILHEIT  AND  BARRE??’ 

: 12]  REFTEXTi-  • WELCH  AND  MIZUSHIMA  ■ 

'13]  RE  EC  EX  7*4*  • EVENSON  AND  MIZUSHIMA  • 

1 1 * J REFTEXTi-' ALBRITTON.  HARROF.  SCHMELT EKOPF , 4 SO  Z4S6” 

115]  • tUJi-O  0,V  7SiVS  JUBiSS  //  4V0  JT.  4WO  00  /4LOS  O.V  P.  118. 

[15]  REFTEXTt-'TOHJTA.  MIZUSHIMA.  HOWARD.  AND  EVENSON' 

[17]  ■ S45£0  O.V  PS//47P  COMMUNICATION  FROM  M.  MIZUSHIMA  TO  M.  GREENEBAUM,  7/7/75. 

[18]  'EL  EASE  TYPE  REFERENCE  NUMBER' 

[19]  X-HEFNO-J 

[20]  JU./ATSUU] 

[21]  Sl*/oTSl[jr] 

[22]  BI-VSCBIlX] 

[23]  L4V0*/£c.T4.V0U  ] 

[24]  04.Vl*/£O,4.VlUf] 

[25]  M’JO-VECMUOlX  ] 

[26]  MJl-VECMJllX] 

[27]  REFTEXTX-i  42o’  • 

[ 28]  REFTEXTXll  ■.  id REFTEXTi  1-REFTSl I?1 
[ 29]  REFTEXTXl  2;  ioS£T7£ir72]*.l?£'P7F*77 
[ 30  ] REFTEXTXl  3 ; i oREFTEXTi  1-REFTEXT3 
[ 31  ] SSfTSJTJ’VC  4;  xoREFTEXTn  l-REFTEXT* 

[32]  REFTEXTXli;  icREFTEXTSl-REFTEXTS 
[ 33  ] P.?f7ar?m;  \oREF7EXT6l-REFTEX76 
[34]  • THANE  YOU. 


[35]  "RE  PARAMETERS  ACCORDING  TO  ' -.REFTEXTXl  X -.  ];  • 
iR'F.  '-.X,')  AREAS  FOLLOWS: 


[36] 

•so  * 

■iSO;1  OSZ;  • 

[37J 

■si  * 

' ;S1  i ' GHZ-.  ' 

[33] 

•S2  = 

i 

•;S2;  • GHZ'. 

[39] 

•04. VO 

- ' iLAMOi  ' GHZ-. 

[40] 

•04.V1 

• 

- ' :LAM\  ; • GHZ-. 

[41  ] 

•vyo  = 

’-.MJO;  ' GHZ-. 

[42] 

•W1  s 

'-.MUU'  GHZ. 

[43] 
[44  ] 

[45] 
[4  6] 


'TEMPERATURE 
' "TEMP-l 
TEMP-  29  6 
• REVISED 


IS  ASSUMED  TO  BE  29  6 K-.IF  THIS  IS  INCORRECT. 
)"  WHERE  DESIRED  NUMBER  IS  INSERTED.' 

10  JULY  1975 


WRITE:  ' 


Fig.  C-l:  Listing  of  APL  function  PARAMETERS 
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V PaRAMSTE  INUPPER  -,X  ; VECBO'.VECBX ; VECB  2 ; VECLAM  Oi  VECLAM! ; VECMU  O'.  VECMUX  ‘.RF.FTEXTS'.VECGV 


[1] 

ft  ' REFNO ' IS  5,  CORRESPONDING  TO  'REFTEXT  5';  ' V ' IS 

0 OR  1 . 

[2] 

ft  CORRESPONDING  TO  V"  . ( ONLY  V'*V''  TRANSITIONS  CONSIDERED 

[3  ] 

VECB 0*-(  1 .437708  1 . 421914  ) *SPEEDOF LIGHT 

[4] 

VECBX*-( “4 . 84£~6  "4 . 840 2 E~ 6) *SPEED0 FLIG H T 

[5] 

VECB2*-  0 0 

[6] 

VECLAM 0*1. S*(l.3239  1 . 32  87  ) *SPEEDOFLIG  HT 

[7] 

VFCLAM 1-0  0 

[8] 

VECMU 0*-( “8  .436  -8 . 9 7 ) *0 . 0 01  *SPE EDO FLIG HT 

[9] 

VECMU  l«-0  0 

[10] 

VECGV-0  1556.378ft  INVERSE  CM 

[11  ] 

REFT  EX  TS+- 'ALBRITTON,  HARROP . SCHMELT  EKOPF , AND  ZARE 

(REF.  5) 

FOR  V 

ii  - y i i i i - i 

[12] 

VV*-X  •, REFNO*-  5 

[13] 

X+VV+X 

114] 

BO-VECBOlX ] 

[15] 

BX-VECBllX  ] 

[16] 

B2*-VECB2lX  ] 

[17  ] 

LAMO-VECLAMOlX  ] 

[18] 

LAMl-VECLAMllX  J 

[19] 

A/0  0«-/i’0Wi/0[lf  ] 

[20] 

A/Ol-yffCWOlU] 

[ 21  ] 

0/«-/£C0/[X  ] 

[22] 

TOff  PARAMETERS  >!/?£:  AS  FOLLOWS  ACCORDING  TO 

' iREFTEXT  5 • : 

1 

[23] 

•BO  - '-.BO-,'  GHZ-,  ' 

[24  ] 

•Bl  = ' iB  1;'  GHZ',  • 

[25] 

•B2  = ' ',B  2;  ' GHZ-, 

1 

[ 26] 

'LAM 0 = '\LAM0-,'  GHZ',  ' 

[27] 

•LAMX  = • ; Li4Wl  ; ' GHZ', 
1 

[28] 

•WOO  = ';MU 0; ' GOZ;  ' 

[ 29  ] 

•W01  = ' ; WOl ; • GHZ. 
1 

[30] 

'10  OOiTS  OF  INVERSE  CM,  THE  ABOVE  PARAMETERS  ARE:' 

[ 31  ] 

+PARAMVEC*-(B0,BX,B2,LAM0 ,LAMX ,MU 0 ,MUX ) *SPEEDOF LIGHT 

[32]  • TEMPERATURE  IS  ASSUMED  TO  BE  296 K',IF  THIS  IS  INCORRECT,  WRITE'.' 

[33]  ' "TEMP*-(  )"  WHERE  DESIRED  NUMBER  IS  INSERTED.' 

[34]  TEMP*-  29  6 

[35]  X«-ln  CALCULATION  FOR  016=016  ONLY,  FOR  V'*V"=X  ONLY. 

[36]  ft  THIS  IS  A MODIFIED  FORM  OF  'PARAMUPPER'  FOR  USE  WITH  'A  BCD  ' OR 

[37]  • 'ABCDLP'  ( THEORY  0F  W.  R.  STEINBACH  INSTEAD  OF  TINKHAM  AND  STRAND  BERT ) 

[38]  » VERSION  OF  28  JULY  1975  . 

V 


Fig.  C-2 : Listing  of  APL  function  PARAMSTEI NUPPER 
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APPENDIX  D 


Temperature  Dependence  of  Rotational 
State  Sum  for  Molecular  Oxygen 


ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=016 


EMPLOYING  THE  PARAMETERS  OF  REF.  1 (V  = 0) 


TEMP 

EXACT  SUM 

77  K 

56.64467 

100K 

73.32311 

125  K 

91 .45851 

1 SOK 

1 09 .59893 

1 7 5/C 

127. 74352 

200K 

145. 89184 

■>7  5 K 

164 .04365 

2 -OK 

1 ^2 .1  9881 

273  K 

198.  90444 

27  5/C 

200.35724 

296  K 

21 5.61279 

300  K 

21 8 . 51 886 

32  5a' 

236 .68364 

CLASSICAL 

RATIO 

5 5.8341  6 

1. 01452 

72.51190 

1.01119 

90. 63  987 

1.00903 

1 08.  76785 

1.00764 

1 26.89  582 

1, 00668 

145.02380 

1.00599 

163.15177 

1.00547 

181. 27975 

1.00507 

197. 95748 

1. 00478 

199.40772 

1 .00476 

214.63522 

1. 00455 

217.  53  570 

1 .00452 

235.66367 

1 . 00433 

ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=016 


EMPLOYING  THE  PARAMETERS  OF  REF . 2 (/  = 0) 


TEMP 

EXACT  SUM 

CLASSICAL 

RATIO 

11 K 

56.64  683 

55.8371 6 

1.01450 

100  K 

73.32569 

72 .51 579 

1.01117 

125  K 

91 .46140 

90 . 644  73 

1 .00901 

1 50  K 

109 . 60200 

108.77368 

1.00762 

11  5K 

1 27.74661 

126.90263 

1.00665 

200  K 

145.89481 

14  5. 031  57 

1. 00595 

225  K 

1 64 . 04635 

163.  160  52 

1 .00543 

2 50  K 

1 82. 20110 

181.28947 

1 . 00503 

213K 

1 98.90622 

1 97. 9681 0 

1 .00474 

27  5K 

20  0. 3 5896 

1 99 .41  841 

1. 00472 

296  K 

215.61392 

21  4.64673 

1 .00451 

300  K 

218.51987 

217 . 54736 

1 . 00447 

325  K 

236 .68380 

23  5.67631 

1 . 00427 

t-1/306-3-14 
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ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN . 016=016 
EMPLOYING  THE  PARAMETERS  OF  REF.  3 (/  = o) 


TEMP 

EXACT  SUM 

77  K 

5 6.64  74  0 

1 OOA 

73.32662 

1 25* 

91 .46282 

1 5 0 A' 

1 09.  60401 

1 7 5 A 

1 27. 74932 

2 OOA 

145.89832 

2 2 5 A 

1 64 .05077 

2 5 OA 

1 82. 20653 

273  A 

1 98.91267 

275A 

200.36551 

296A 

21 5. 62149 

3 OOA 

21 8 . 52765 

3 2 5 A 

236.69291 

CLASSICAL 

RATIO 

55.83725 

1 . 01451 

72  . 51  591 

1. 01118 

90 . 64488 

1 .00902 

108.77386 

1.00763 

126. 90284 

1.00667 

145.03181 

1.  00  597 

163. 16079 

1.00545 

181.28976 

1. 00506 

197. 96842 

1 .00477 

199.41 874 

1 . 00475 

214 . 647  08 

1 .00454 

217 . 54772 

1.00450 

23  5.67669 

1 . 00431 

ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN.  016=016 
EMPLOYING  THE  PARAMETERS  OF  REF.  4 (V  = 0) 


TEMP 

EXACT  SUM 

77  A 

56.64740 

1 OOA 

73 . 32662 

12  5A 

91 . 46282 

1 50  A 

1 09.60401 

175A 

127.  74932 

2 OOA 

145. 89832 

2 2 5 A 

164 .05077 

250  A 

182.20653 

2',3  A 

1 98.91  2 68 

275A 

200.36551 

29  6 A 

21 5. 62150 

3 00  A 

213.52765 

32  5A 

236. 69291 

CLA SSICAL 

RATIO 

55.83725 

1. 01451 

72 . 51591 

1 .01118 

90.  64488 

1.00902 

108 .77386 

1 .00763 

1 26. 9 0284 

1 . 00667 

145.031 81 

1.00597 

1 63 . 1 6079 

1.  00545 

1 81. 28975 

1. 00506 

1 97.  96842 

1. 00477 

199.41 874 

1 . 00475 

214.64708 

1.00454 

217. 54772 

1 . 00450 

235.67669 

1. 00431 
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ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=016 
EMPLOYING  THE  PARAMETERS  OF  REF.  5 (V  = 0 ) 


TEMP 

EXACT  SUM 

CLASSICAL 

RATIO 

HE 

56. 64484 

5 5 . 83  61  0 

1.01448 

100  E 

73 . 32372 

72.51 442 

1.01116 

12  5 E 

91.45955 

90.64302 

1 .00901 

1 50  E 

1 09. 60038 

108.77163 

1 .00762 

175K 

1 27.  74  532 

1 26 . 9 0023 

1 . 00666 

200  E 

145. 85396 

14  5.0  2884 

1 .00597 

225  E 

164.04605 

163.1 5744 

1 . 00  54  5 

250  E 

182.20144 

1 81 . 28605 

1 . 0050  5 

21 2 E 

1 98 . 90725 

1 97. 96436 

1 .00476 

27  SE 

200.36006 

1 99.41  465 

1 .00474 

296  E 

215.61 574 

214.64268 

1 . 00453 

300  E 

218.52183 

217  . 5432  6 

1 . 004  50 

32  5 E 

236. 68672 

235.67186 

1 . 00431 

ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=016 
EMPLOYING  THE  PARAMETERS  OF  REF.  6(^=0) 


TEMP 

EXACT  SUM 

CLASSICAL 

RATIO 

HE 

56 . 64749 

55.83732 

1 . 014  51 

1 00  E 

73.32675 

72 . 51 600 

1. 01118 

125  E 

91  .4  6299 

90.64  500 

1 .00902 

150  E 

1 09.6  0423 

108.77400 

1 . 00763 

175  E 

1 27.749  58 

126.90300 

1 .00667 

200  E 

145.89863 

^145.03200 

1 . 00598 

22  5AT 

1 64 . 0 51  14 

163.16100 

1 .00546 

250  E 

1 82. 2069  5 

181  . 29001 

1 . 00506 

213E 

198.91315 

197.96869 

1 .00477 

212E 

200.36599 

199.41 901 

1 .00475 

296  E 

21 5.62202 

21  4.64737 

1 . 00454 

300  E 

218.52818 

217 . 54801 

1 . 00451 

325  E 

236.69350 

235.67701 

1 . 00431 
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ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OX  KEN,  016=016 
employing  the  parametehs  OF  REF.  7 ( v = 0) 


TEMP 

EXACT  SUM 

11 K 

56 . 64  74  8 

i oo  a: 

73.32672 

12  5# 

91 . 4 629  5 

1 50# 

109 .6041 6 

17  5# 

1 27.  74949 

200# 

145.89851 

2 2 5# 

164.05098 

2 5 0# 

182.20676 

273# 

1 98. 91292 

275# 

200.36575 

2 96# 

21  5. 62175 

3 0 0# 

21  8 . 52791 

32  5# 

236.69317 

CLASSICAL 

RATIO 

55  . 83732 

1 . 014  51 

72.  1600 

1 . 01 1 1 8 

90.64500 

1.00902 

1 08.774  01 

1 . 00763 

126.90301 

1.00667 

14  5. 03201 

1 . 00597 

163.16101 

1 . 00  54  5 

181. 29001 

1.00506 

197. 96869 

1 . 00477 

1 99 .41  901 

1 . 00475 

214. 64737 

1 . 004  54 

217  . 548  01 

1 . 004  50 

235.67701 

1 .00431 

ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=018 
EMPLOYING  THE  PARAMETERS  OF  REF.  7 (V  = 0) 


TEMP 

EXACT  SUM 

77# 

.11  6 . 284  88 

1 00# 

1 51 . 6071 2 

125# 

190.01360 

1 50# 

228.42996 

175# 

266 . 854  66 

200# 

305. 28693 

225# 

343. 72634 

250# 

382.17263 

273# 

417. 5491 5 

27  5# 

4 20 . 62564 

296# 

4 52.  931  28 

300# 

459.08526 

32  5# 

4 97.  55143 

CLASSICAL 

RATIO 

118.23962 

0. 98347 

1 53  . 5 579  5 

0 . 9873  0 

1 91  . 94744 

0 .98993 

230. 33693 

0 . 991  72 

26  8 . 72  64  2 

0 . 99303 

307.11590 

0.99404 

34  5 . 50  539 

0. 99485 

383 . 89488 

0.99551 

419.21321 

0.99603 

422 . 2 84  37 

0.99607 

4 54 . 531  54 

0 .99648 

4 60 . 67386 

0. 99655 

499.06334 

0.99697 
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ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN , 018=018 
EMPLOYING  THE  PARAMETERS  OF  REF.  7(^=0) 


TEMP 

E*>icr  SOM 

CLASSICAL 

RATIO 

11 K 

63. 28076 

62  .81336 

1 . 00744 

1 00* 

82.044  95 

81  . 5 75  80 

1 . 00  575 

125  K 

102.44779 

101.96975 

1 . 00469 

1 50* 

122.85599 

12  2.3  637  0 

1 . 00402 

175* 

1 43 . 26871 

142.75765 

1 . 003  58 

200* 

1 63 . 68549 

163.15160 

1 .00327 

225  * 

1 84. 1 061 2 

1 83  . 54  5 5 5 

1 . 00305 

250  * 

204. 53043 

203.93950 

1 . 00290 

273* 

223.32398 

222.70193 

1 . 00279 

2 7 5* 

224.95834 

224.33344 

1 . 00279 

296  * 

24  2.120  53 

241 .4  64  36 

1 . 00272 

300* 

245.38980 

244 . 72  739 

1 .00271 

32  5* 

265.82476 

265.12134 

1.00265 

ROTATIONAL  PARTITION  FUNCTION  FOR  MOLECULAR  OXYGEN,  016=016 
EMPLOYING  THE  PARAMETERS  OF  REF.  5 ( V = 1 , 


TEMP  EXACT  SUM  CLASSICAL  RATIO 


77* 

57. 22604 

56 .4  5631 

1 . 01363 

100* 

74 .09050 

73. 31  988 

1 .01051 

12  5* 

92.4  281  5 

91 . 64  985 

1 . 0084  9 

1 50* 

110.77088 

1 09  < 97982 

1 .00719 

175* 

1 29.11  783 

1 28.30979 

1 . 00630 

200* 

1 47.4  68  58 

146.63  976 

1 . 00  565 

2 25  * 

165.82288 

1 64 . 9 69  73 

1 . 00517 

250* 

1 84.1  80  59 

183. 29969 

1 . 00481 

273* 

201 .07261 

200. 1 6327 

1 . 004  54 

275* 

202. 541 61 

201. 62966 

1 .00452 

296* 

21  7.  96740 

217. 02684 

1 . 00433 

300* 

220.90590 

21  9 . 959  63 

1 .00430 

32  5 * 

239.27341 

238. 28960 

1 .00413 
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appendix  e 


Line  Strengths 

Computed  hy  ,,s  at  ^ and  gQfiK  with  -s , 

Results  Given  in  Ref.  8 

text,  pp.  27  and  34  for  discussion  of  this  Appendix, 
was  computed  using  the  APL  program  GBBSTRENGTH.  It  is 
assumed  that  the  cm-atm  quoted  in  Ref.  8 refers  to  273K  and 

1 atm,  so  that  1 (cm-atm)^  = 2.686754(84)  El9  molecule  cm’* 
Loschmidt's  constant,  L o). 13,23 

Th*  *««  L0  used  in  the  APL  program 
GBBSTRENGTH  should  have  been  replaced  by  Lo/0. 99519  in  computing 
the  integrated  line  strength  in  units  of  cm-  per  (cm-atml 

nee  p.  17  of  text).  The  0.5*  change  in  our  computed  line 
strengths  does  not  affect  the  conclusions  to  be  drawn  from 

the  comparison,  i P *-v,  = +-  *-u 

*,  hat  the  results  of  Ref.  8 may  be  in 

error  by  as  much  as  a factor  of  2.  Note  also  that  the  APL 
program  GBBSTRENGTH  was  used  only  in  generating  this  Appendix- 
the  line  strengths  elsewhere  in  this  report,  given  in  units  of 

nm  ^ 1 - — *■» 


cm 


per  molecule  cm  2,  have  been  calculated  correctly. 
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TE-P  • 273 KJ  <T  • 119*7  INVERSE  CM)  X ■ 1 

REP’.O  • 7 

N1N»UT!  1 3 5 7 9 11  13  IS  17  19  21  23  15  27  29  31  33  35  37 

39  *1  93  95  97  99  51  53 


the  line&trenqths  op  the  submillimeter 

units  CP  l*OE>6  INVERSE  CM  PER  CH-AT* 


TRANSITIONS  OP  016-01** 
AT  STP>«  ARE  RS  POLLOnS 


IN 


sp: 


0 • 699*  1*959  2*S03  3*152  3*057  2**S 

0 • J66 1 0*1981  0*1002  0*09799  0*02105 

0**009336  0*000190*  9.293E-5  1*233E-S 

***22E-S  9*86SE-9 


2*09  1*519  1*015  0*631* 

0*0087*3  0*00392*  0*001257 

3*338E>*  8 *5E-7  2*09lE-7 


so: 


7*93*  12*37  15*19 

1*597  0*858*  0 • *32 

0*0005931  0*0001807 

* • 12E-B 


15*79  1*«*B  12*33  9*52S  6*801  **509  2*777 

0*203*  0*09  0*03733  0«01*5*  0*00532*  0*001832 

5* 181E-S  1.399E-5  3*S*IE-6  B*S9lE-7  1*9311-7 


SMi 


3*2*  9*261  9*75*  9*71  **228 

0*9269  0.2279  0*1191  0*05357 

0*0009752  0*0001535  9***8E-S 

**95*E-8  1 • 0S6E-8 


3**79  2*493  1*8*2  1*221 

0*02359  0 • 009757  0*00379 

1 • 337E-5  3 • 60RE-*  9*159E-7 


0*7965 

0*001389 

2.199E-7 


QEBBIE*  BURROUGHS*  AND  BIRO*  PROC*  ROT  • SOC*  I LONDON 


579 

O 

m 

o 

*- 

119*9)  GIVE  POR 

sp: 

0*62 

1*5 

1*8  1*8 

so: 

5.7 

8*7 

9*5  8*3 

SH 

2*1 

2*9 

2*9  2** 

THE  PIRST  SEVERAL  LINES  I 

1** 

*•3 

1*8 


A*  VOL*  310*  PP* 

N THE  SAME  UNITS) I 


RATIOS  OP  OUR  RESULTS  POR  K • 1 THAQUOP  9 TO  ThOBE  OP  OCBBIE*  ET  AL I 


1*128 

1*30* 

1 *557 

1*751 

2*188 

1*3* 

1**22 

1*593 

1*902 

2*325 

1*553 

1«*69 

1*639 

1*962 

2*3*9 

T-l/3 06-3-14 
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TE*'P  • 2yb«)  XT 

RFTNO  • 7 

■ 

205.7  inverse  oj 

X ■ 1 

Mnput:  1 3 5 

7 

9 11  13  15  17 

19  21  23 

26 

27  29 

31  33  35  37 

J 9 *1  -3 

-5 

-7  -9  51  53 

Tnfc  - INE-jTRENGThS 

OF 

The  SubMlLLIMETER 

TRANSITIONS 

OP 

016*016/ 

IN 

-MTS  OF  l*0E-6  INVERSE  CM  P£R  Cm-AT-  IaT  STP)/  ARE  AS  FOllOwS! 


SF:  C.-JS77  1 *687  2 ••mb  2*799  2*775  2*-7l  2*011  1*511  1*055  0*6B6I; 

0 . **  1 85  0*2391  0*  1283  0*06-77  0*03078  0*01878  0*00582  0*0023l9 

0 * w0087£5  0*0003101  0*00010-1  3*307E-5  9*983E>6  2*B23E>6  7.596E-7 

1 * 935E-7  - • 6 7 1 £ ■ 8 

so:  6*321  10*65  13*2  1-.01  13*29  11*5  9*165  6*781  -*678  3*018  1*825 

1*036  0*553  0*2778  0*1315  0*05869  0*02-71  0*009818  0*003685 

0 • uO 1 307  0*000-38  0*0001389  -*165E>5  1*1821-5  3.176E-6  8*083E-7 

1 *a-9E-7 

Sh:  2*786  3*67  -*l-8  -*182  3*838  3*2--  2*563  1*857  1*268  0*8115 

0.-87-  0.2751  0*1-61  0*07312  0*03--9  0*0183-  0*006--  0*002553 

0 • u009559  0.0003383  0*0001132  3*S83E»S  1.073E-5  3*0-lE-6  8*16E>7 

2 » u 7-E-7  - • 995E- 8 


QEBBIE/  oURRCuOmS/  AND  BIRD/  PROC*  ROT*  SOC*  I LONDON  I A/  VOL*  310/  PP* 
579  TO  590  <19691  OlVE  F 0*  Th£  FIRST  SEVERAL  LINES  ( IN  THE  SAME  UNITS)! 


sf  : 

0*62 

1*5 

1*8 

1*8 

1 *- 

so: 

5.7 

8*7 

9*5 

8*3 

6*3 

sm  : 

2*1 

2*9 

2*9 

2*- 

1*8 

RATIOS  ON  OuR  RESulTS  FOR  < ■ 1 TmRQuQm  9 TO  TnOSE  OP  QEBBIE/  ET  AL i 


f : 

0**6- 

1*125 

1*359 

1*555 

1*982 

g: 

1*1-- 

1*226 

1*39 

1*688 

2*11 

h : 

1*326 

1*266 

1.-3 

1*7-2 

2*132 

i 

I 
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APPENDIX  F 


RRI  Absorption  Line  Parameters 
for  Molecular  Oxygen  Isotopes 
J-60160.  l60l80t  and  l80le0  ( AFCRL  Format) 

In  the  following  eight  file  listings,  the  contents  are 
labelled  at  the  top  of  each  page  with  the  file  names  as  well 
as  brief  descriptions  of  the  contents.  The  units  are  as  follows 

FREQ:  Transition  frequency:  (cm-1) 

STRENGTH:  Line  strength  at  296K:  (cm  1/molecule  cm-2) 

WIDTH:  Line  half-width  at  half-maximum:  (cm-1  atm-1) 

E'':  Energy  of  lower  state  of  the  tran- 
sition, with  respect  to  the  ground 
state,  including  AG1/2:  (cm-1) 

The  sequence  of  the  eight  file  listings  is  as  follows: 
AF76OKI5VO:  Microwave  fine  structure  lines  of  l60l60,  v = 0. 
AF7KI5VO:  Submillimeter  rotational  lines  of  l60l60,  v = 0. 
AF76OK25VO:  Microwave  fine  structure  lines  of  l60le0,  v = 0. 
AF7K25VO:  Submillimeter  rotational  lines  of  l60l80,  v = 1. 
AF76OK35VO:  Microwave  fine  structure  lines  of  le0l80,  v * 0. 
AF7K35VO:  Submillimeter  rotational  lines  of  l80le0,  v = 0. 
AF56OKI5VI:  Microwave  fine  structure  lines  of  160l60,  v = 1. 
AF5KI5VI:  Submillimeter  rotational  lines  of  l60l60,  v = 1. 
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«M  ABSORPTION  1 2 n £ P6PA-CTEPS 
*A7io«'.3*0  60,  aii  £T£  • l Pi . . 7, 


AON  MOLECULAR  OXVNEN  ISOTCPtS  IaFC*l  rO»KAT) 
01*»01i  *ICNO**vE  LINES,  VNu»EME  *IDTnS,  v • s 


2 

3 


* 

7 

a 

9 

:o 

:i 

12 

13 

1* 

19 

1* 

1 7 

ia 

19 

20 
21 

22 

23 

2** 

23 

2* 

27 

2a 

29 

30 

31 

32 

33 
3* 
39 
3a 
37 
3a 
39 
* 0 
*1 
*2 
*3 
*6 
*3 
*1 

*7 


‘"to  5VP£NQT*i  »]OTh  £1*  y 1 J|  n 1 


1 *5*117 

2*  ol£>33 

*032 

1 *51*99 

8*53£-33 

• 032 

1 • 511*3 

3 * 36£*32 

•032 

1*39171 

1 *23E*31 

•032 

1 *91911 

* * * IE-3 1 

•032 

1*13235 

1 • • 7£  <•  30 

• 032 

1*1*933 

* • 12E*30 

• 032 

1 •11135 

1*38E*29 

•032 

1 *11362 

3*S7£«29 

•032 

1 • 7Sl7i 

1 *03E>28 

•032 

1 • 71791 

2*9fl£-2« 

• 032 

1*7332* 

1 • 09E-2N 

•032 

1 • '5212 

1 • 31E-2  7 

•032 

1*77012 

2 • 83E-2  7 

•032 

1*71776 

3 • 13E-2  7 

•032 

1 *80531 

1 *09E-21 

• 038 

1 *82363 

1 *83E-26 

•039 

1 *1*199 

3*CC£«21 

•037 

1*81073 

**60E-21 

• 038 

1*87675 

2*76£-26 

•0*3 

1 *81001 

6 • S8Ea21 

• 038 

1 *90026 

8*77£«26 

• 039 

1*92175 

1 *08E>23 

•0*1 

1 *9*3*8 

1 *22E-23 

• 0*3 

1*9*537 

7*55E-21 

• 0** 

1*97331 

l*2lE-23 

• 0*6 

1*9877* 

1 • 1 1 £*29 

•0*2 

2*01159 

1*1  3£*23 

• C** 

2*01919 

1 • 31E-25 

• 0*1 

2*03976 

1 *33E*2B 

• 0*0 

2*061*3 

1 *231*23 

•039 

2*08181 

1 • 03E -25 

•038 

2*08*32 

8**1£>21 

• 0*7 

2*10139 

8 *23E-21 

• 03* 

2*120*2 

3*97E*21 

• 031 

2*13907 

* • 0*E>26 

•039 

2*137*6 

2 • 31£*21 

• 039 

2*1736* 

1 *32E*21 

• 032 

2*19318 

8 • *9E*27 

•032 

2*21110 

***3E*27 

•032 

2*229*3 

2*20E*27 

•032 

2*2*720 

1 • 02E*2  7 

• 032 

2*26*92 

* • *8E*28 

•032 

2*28210 

1 • 8SE*2S 

• 032 

2*30021 

7 * ?*E  *29 

• 032 

2*31789 

2 • 17E-29 

• 032 

2*33331 

9 • 29E-30 

• 032 

*07*. 332 

0 

93 

98 

3778*223 

0 

91 

31 

3*92*781 

0 

*9 

*9 

3218*270 

0 

*7 

*7 

293*. 737 

0 

*3 

*8 

2702.223 

0 

*3 

*3 

2*10.77* 

0 

*1 

*1 

2230. *23 

3 

39 

39 

2011.213 

0 
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RIVERSIDE  RESEARCH  INSTITUTE 


A8SCRRTI0N  LlN£  »A«4f£T£HS  f 0 " 
*'340<!5vl  80/  81/  ETC.:  R£f.  5,  ci*-Olfc 


"QlECulA*  OXTQEN  ItOTOR£g  (AfCRL 

* I CRO** VC  LINES/  KRuRENIE  »IOThS/ 


'"CO  STR£nQTh 

*•  I'lllU  1 • 6 3E • 3 3 

*9  2.399*3  5 • 83£ • 3* 

30  ?.*1743  1«75£*3* 

31  2*43585  *.9*E«33 

32  2 • *3*02  1 • 32E -35 

33  2 • * 72 1 • 3 • 3AE-34 

3*  3.97713  3 • 1 7£  -29 


" IOT«  £'*  V*  j'X* 

•032  *221.0*5  1 *3  *3 

•032  **77.73*  1 ,5  *S 

•032  *738.297  1 */  *7 

•032  3C09 . 7 1 1 1 „*  *9 

•032  5291*932  1 51  3l 

•032  55|4 . 9 1 1 1 53  53 

•050  153*  * 3 78  1 1 i 


v"  j>  >«>  > jo  date 

**  *3  *3*  73 

*4  43  *3*  73 

48  *7  *7*  73 

SO  *9  *9*  73 

32  31  31*  73 

34  33  33*  73 

Oil-  73 


rORNAT) 
V • 1 
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44  7 

44  7 

44  7 

44  7 

44  7 

44  7 

44  7 
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»•*!  *(!Swfc«f,T  ICS  lI*»E  PARAMETERS  FOR  HC„ECk,LAR  OXTOEN  ISOTOPES  < *FCRt  FQB"AT I 


AF**!*: 


1 

2 

3 

* 

3 

6 

7 

I 

9 

1 0 
n 
12 
13 
1* 
13 
16 
17 
II 

19 

20 
21 
22 
23 
26 

25 

26 
27 
26 

29 

30 

31 

32 

33 
36 

30 

36 

37 

31 
39 
60 
*1 
62 
63 
66 
•5 
66 
67 


bCj  61*  ETC'S  REF.  5*  016-016  SUB***  i-lsEB#  xRuREME  »1DT*6*  v • 1 


FREu 

STR£nGT* 

« JOT*. 

l2*13ilB 

1 . 1 *E-29 

* 0*8 

16.02102 

1 .25E-28 

• 0*5 

1 6 > 10831 

5 * 36E-29 

•0*5 

23*57570 

3 • 2-E-29 

• 0*0 

£5*5377* 

2 • 05E-28 

<0*. 

£7.55283 

7 • 08E-29 

*066 

J6. 98265 

* * 7 1 E «29 

• 0*3 

J6  * 9 B 309 

2**5E-28 

• 0*2 

36 « 95958 

8 • 0 1 E-29 

• O', 2 

66*37361 

5 • -0E-29 

• 0*2 

.8.60290 

2 • 7 1 E-28 

• 0*1 

30*35056 

8.09E-29 

• 0*1 

57.75231 

5.37E-29 

• 0*1 

39*80650 

2 • 58E  *28 

• 0*0 

bl  *7296. 

7 • *5E-29 

• 0*0 

t,9.11929 

* • 80E-29 

•0*1 

71  * 19601 

2 • 2*E -28 

*039 

73*09662 

6 • 32E  *29 

• 039 

«0. 67360 

3 * 92E -29 

• C 38 

82.57138 

1 • 79E-28 

•038 

B 6 * 65053 

* • 9 7E  *29 

• 038 

9K61322 

2.96E-29 

• 036 

93*93166 

1 -33E-28 

• 03* 

95*79035 

3 • 65E  *29 

• 03* 

lu3* 13716 

2 • 06E*29 

• 037 

105*27552 

9 • 23E-29 

*036 

107*11629 

2 • 50E-29 

• 036 

1 16.66353 

1 • 36E-29 

*036 

116.60162 

5 • 99E-29 

*035 

1 1 8 * 62066 

1 *61E*29 

• 035 

125*73056 

8*35E*30 

• 035 

127.90771 

3 • 6*E-29 

• 035 

129*70769 

9*  73E-30 

*035 

1 36*99667 

6 . B0E-30 

• 035 

139*19266 

2 • 08C-29 

*032 

160*97360 

5*53E*30 

*032 

168*23965 

2*602-30 

• 032 

150*65651 

1 • 12E-29 

• 032 

152*21656 

2 • 96E*30 

• 032 

159*65767 

1 • 32E-30 

• 032 

161*69168 

5.68E-30 

*032 

1 b3  * 63680 

1 .69E-30 

• 032 

170*66931 

6*361-31 

*032 

1 72*90175 

2*7l£*30 

• 032 

1 76.626*6 

7 • 1 1C  — 31 

• 032 

181*81251 

2 • B6E*  31 

• 032 

166*08351 

1 *22E*30 

• 032 

E ' 1 

v 1 

J' 

X ' 
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1 

2 

3 

1558*665 

1 

( 

3 

1558**65 

1 

3 

a 
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1 

6 

8 
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6 

9 
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5 
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10 
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1 
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18 
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17 
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1 

16 

17 
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17 
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1 

18 

19 
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18 
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21 
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21 
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21 

21 
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«USOrfPT'ON  LJN£  8*8*H£T£HB  T Q*  HQlECuL**  OXTQEN  I80T08ES  I 6ECNL  EO*H*T> 


*EJ« 1 5v l O0< 


E«EU 


*8 

U5 

78960 

*9 

192 

9ojoo 

50 

195 

23o93 

51 

196 

92257 

52 

2jo 

0o62b 

53 

2w6 

35ol7 

5* 
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02338 

55 
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5* 
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03938 

57 

?1  9 

09022 

58 

?26 

1 oo  1 1 

59 
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*0 

2J0 
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61 
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62 

2 j9 
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63 

2o  1 

1 1058 

60 
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09127 

65 
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66 
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67 
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68 
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<*031o 

69 

26? 
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70 
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71 

2/2 

29056 

72 

2/3 

8500O 

73 

2«0 
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7 0 

2n3 

13286 

75 

2»o 

67o  1 o 

76 
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0/017 

77 

293 

92619 

7a 

295 

05130 

79 

3s#2 

20252 

80 

3jo 

6 7o  70 

81 

3w6 

17965 

B 1 # ETC.l  N£E. 


STNfN3T«  »;otm 

3 • 1 9£«  3 J <032 

1 • 22£»  j 1 .032 
5 • 1 9t * 31  *032 
1.35E-31  .032 
o*92£«32  .032 
2*08£«31  .032 
5**2£«32  >032 
1-B7E-32  .032 
7 • 91 £“32  .032 
2.C3E-32  .032 
6*J?E-33  .032 
2 « 8o£«32  .032 
7-35E-33  .032 

2 *29E«33  .032 
9.63E-33  .032 
2**9E-33  .032 
7 • 36E“3o  .032 
3*09E-33  .032 
7.9BE-3*  *032 
2*2*E»3o  *032 
9 « oO£»3o  .032 
2**>2E-3<*  -032 
6*06E-35  .032 
2*71£«30  -032 
6 « 96E*35  .032 
1*76E-3B  .032 
7.3BE-3S  .032 
1.90E-3S  .032 
* • 56E  >36  «03i. 
1-91E-35  .032 
6*89£-36  -032 
1*12E»36  .032 
o*67£-36  .032 
1 • 20E»  36  .032 


* 01 6*0 1 6 SuBn" 


£ ' ' v • j • 
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1 

33 
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1 

30 

iloo.516 

1 

30 
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1 

35 
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1 

36 
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1 

36 
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1 

37 
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1 

38 
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1 

38 
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1 

39 
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1 

oO 
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1 

00 
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1 

ol 
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1 

o2 
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o2 
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1 

o 3 
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1 

*,o 
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* 

0* 
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1 

o5 

OO80* 1 3o 
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o 6 
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1 
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1 
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08 
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o9 
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1 
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51 
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1 

52 
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1 

5? 
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1 

51 
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1 

50 

558*  * 9 1 1 

1 

5* 

558* • 9 1 1 

1 

55 
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APPENDIX  G 


for  Molecular  Oxygen  Isotopes 
160160.  l6olao.  and  l80le0  Who s e 


See  text,  page  33  for  discussion  of  the  relationship  of 
the  file  OXYGENEXIST  to  the  eight  files  listed  in  Appendix  F. 
The  units  are  the  same  as  in  Appendix  F. 


RIVERSIDE  RESEARCH  INSTITUTE 


J,,0E''“ST  1!OIO,ES  ,.,„E 

A < f . 


1 

2 

3 

4 

3 

6 

7 

8 
9 

10 
11 
12 
13 
1 * 
13 
16 
17 
IS 

19 

20 
21 
22 
?3 
2* 
23 
26 

27 

28 

29 

30 

31 

32 

33 

36 

33 

36 

37 

38 

39 
*0 
*1 
62 
63 
66 


i!to  -iSsE  !!RKSr2J4EKCEfo03a7f!M  ?IM01i  "IT*  » * » OR  1 

AT  296*.  lInEwIOTiS  I\TE8R0LATEd%b0M  r ?LECUuE  P£R  c*  SO 

“*  s“  •"  I”", 


f«i.4  STRENQTii  -JOT*  £i 


V • 


1 *66953 

1*6661:5 

1 *68362 

1*70074 

1*71796 

1 *73326 

1*73262 

1 • 76631 

1*77012 

1*77254 

1 • 78086 

1 • 76774 

1*7851* 

1*7976« 

1 • 79994 
1*80354 
1*80386 
1*81628 
1*81923 
1*82273 
1*82363 
1*83127 
1 *83877 
1*83946 

1*86199 

1 *86432 
1*83727 
1*85869 
1*86075 
1*86611 
1*87509 
1*87679 
1*87915 
1*88008 
1*88620 
1*88980 
1*83206 
1*89330 
1*90024 
1 *90061 
1 *90302 
1*91282 
1*92175 
1*92293 


<**62E«30 
1 • 3BE-29 
3*87£-29 
1 *03E-28 
2*58E-2B 
6*09£-2B 
1 • 36E-2  7 
3*73E-30 
2*85E-27 
3O3E-30 
'•-9E-30 
3*63E-27 
: ‘C6E-23 
1 *62£«29 
3*5  7£“30 
1 *05E-26 
1*9 1 E-29 
2*56£-29 
9 1 67E-30 
J*32E**29 
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